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We improve and generalize a resummation method of post-Newtonian multipolar waveforms from 
circular (nonspinning) compact binaries introduced in Refs. [1, 2]. One of the characteristic features 
of this resummation method is to replace the usual additive decomposition of the standard post- 
Newtonian approach by a multiplicative decomposition of the complex multipolar waveform him 
into several (physically motivated) factors: (i) the "Newtonian" waveform, (ii) a relativistic cor- 
rection coming from an "effective source", (iii) leading-order tail effects linked to propagation on a 
Schwarzschild background, (iv) a residual tail dephasing, and (v) residual relativistic amplitude cor- 
rections fern- We explore here a new route for resumming firn based on replacing it by its ^-th root: 

l/£ 

Pim ~ firn ■ ^he extreme-mass-ratio case, this resummation procedure results in a much better 
agreement between analytical and numerical waveforms than when using standard post-Newtonian 
approximants. We then show that our best approximants behave in a robust and continuous man- 
ner as we "deform" them by increasing the symmetric mass ratio u = mim2/(mi + 7712)'^ from 
(extreme-mass-ratio case) to 1/4 (equal-mass case). The present paper also completes our knowledge 
of the first post-Newtonian corrections to multipole moments by computing ready-to-use explicit 
expressions for the first post-Newtonian contributions to the odd-parity (current) multipoles. 

PACS numbers: 04.25.Nx, 04.30.-w, 04.30.Db 



I. INTRODUCTION 

One of the prime targets for the eurrently operating 
network of laser interferometer gravitational wave (GW) 
detectors is the inspiral and merger of binary black hole 
systems. To detect and interpret the GW signals from 
such systems one will need accurate templates to match 
theoretically computed signals to the noisy output of 
the detectors. The prime analytical framework allowing 
one to compute (within General Relativity) the GW sig- 
nal emitted by a comparable mass binary system in the 
mildly relativistic regime^ x ^ {v/c)'^ ~ GM/{c^R) <^ 1 
is the post-Newtonian (PN) approximation scheme (see 
Ref. [3] for a review). This raises the issue of the con- 
vergence of the PN expansion, or, in practical terms, of 
the largest value of the PN-expansion parameter x for 
which the currently known PN expansions yield accu- 
rate enough GW templates. Note that, when speak- 
ing of "convergence" in this paper we shall not have 
in mind the mathematical question of whether the full 
PN expansion of, say, the (Newton-normalized) GW 
radiation flux, FTayior^^) ^ J2t=o fni'^'^ogx)x'' is a 
mathematically point-wise convergent series (for some 



^ Our notation is: M = m.i + m2, fJ, = mim2/M, v = li/M = 
m\m,2/{m\ + m2)'^ , O = orbital frequency, v = (GMfl)^^'^, x = 
v^/c^ = {GMQ./c^)'^/^ . We shall generally use x as PN ordering 
parameter, and often use (without warning) units where either 
c = 1 or G = 1. We recall that a term x'" ~ ?;^"/c^" is said to 
belong to the n-PN approximation. 



fixed X belonging to some range) as 71 ^ +00, but 
the more practical question of how small is the numer- 
ical difference (say in the supremum, Leo, norm) be- 
tween the currently known truncated PN expansions, say, 
F^^-yi" = Taylor^ /'(x) = Y.Lg fn{i^^ogx)x^ , for TV = 
3 (3PN approximation), and the "exact" flux F^''^'^'(a;) 
in some physically relevant interval Q < x < Xmax, where 
a^max is equal or close to the value corresponding to the 
Last Stable (circular) Orbit (LSO). We shall then con- 
sider that some resummation method, which transforms 



Taylor , 
TV I 



into FRcsummed^^) (g^^y foj. ^ 



fective if sup^^^. 
icantly smaller than sup^^^ 



^Resummed^^^ ^Exact^^j 



f;^"^'°"(x) 



when Xmax corresponds to the LSO (i.e. 
the extreme- mass-ratio limit ly 0). 



3) is, ef- 
is signif- 

^Exact(^) 

X = 1/6 in 



It was pointed out by Cutler et al. [4] and Poisson [5] 
that the convergence (in the sense just explained) of 
the PN series is rather poor, especially near the LSO 
(i.e. when x ~ 1/6) in the extreme-mass-ratio case that 
they considered. It was then suggested by Damour, Iyer, 
Sathyaprakash [6], to use resummation methods to ex- 
tend the numerical validity of the PN expansions (at 
least) up to the LSO. They used several resummation 
techniques, and in particular Pade approximants. New 
resummation methods, aimed at extending the validity of 
suitably resummed PN results beyond the LSO, and up 
to the merger, were later introduced in the "Effective- 
One-Body" (EOB) approach and used to estimate the 
complete GW signal emitted by inspiralling, plunging. 
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merging and ringing binary black hole systems [7, 8]. 
The BOB method has been recently improved, notably 
by the introduction of a new, resummed, 3+^PN accu- 
rate^ waveform for the £ = m = 2 case [1, 2]. For sev- 
eral, comparable-mass cases, such a waveform (married 
to the EOB dynamics) has been shown to agree remark- 
ably well, both in phase and in modulus, with Numerical 
Relativity data (sec [9] for review of binary black hole 
numerical simulations). For instance, Ref. [10] found a 
phase difference smaller than ±0.025 radians with Jena 
data all over the inspiral and plunge up to merger, while 
Ref. [2] found a remarkable amplitude agreement with 
published Caltech-Cornell data over the inspiral and part 
of the plunge. Let us note in this respect that the use 
of a theoretically less accurate waveform (Newtonian- 
accurate multipolar waveform) still allows for a rather 
small phase difference, but leads to significantly larger 
differences in the modulus [11]. 

The main aim of this paper is to further improve 
the type of resummed multipolar waveform introduced 
in [1, 2] for circularized (nonspinning) compact binaries. 
More precisely, we shall achieve here two goals: (i) on the 
one hand, we shall generalize the resummed £ — m = 2 
waveform of [1, 2] to higher multipoles by using the most 
accurate currently known PN-cxpanded results [12-14] as 
well as the higher PN terms which arc known in the test- 
mass limit [15, 16], and (ii), on the other hand, we shall 
introduce a new resummation procedure which consists 
in considering a new theoretical quantity, denoted below 
as pi.m{x), which enters the (£,to) waveform (together 
with other building blocks, see below) only through its 
£-th power: hira oc {pem{x)Y ■ In this paper we shall pri- 
marily use the small- mass-ratio limit (i^ ^ 0), in which 
one knows both high PN expansions of pem{x) [15, 16] 
and the "exact" value of pimix) from numerical studies 
of test particles around black holes [4, 5, 17], to study 
the quality of the convergence of pI^^°'^{x). Then we 
shall explore the robustness and consistency of our new 
approximants in the comparable-mass case. 

Though we leave to later sections the precise defini- 
tion of the various building blocks of our new, resummed 
waveform, let us already sketch here its structure. The 
basic idea is to write the [i, m) multipolar waveform 
emitted by a circular"^ nonspinning compact binary as 



^ The notation 3+^PN refers to a "hybrid" expression which in- 
corporates both the comparable-mass {v ^ 0) 3PN terms and 
the extreme-mass ratio {y = 0) 4PN and 5PN terms. See below 
for the precise definition of the "hybridization" procedure we use 
here. 

^ In this paper, we focus on the waveform emitted by exactly cir- 
cular orbits (see, however, footnote 12 below). We leave to future 
work the study of "non-quasi-circular" corrections that must be 
introduced in the realistic case of inspiralling and plunging or- 
bits (such corrections have already been introduced in the EOB 
approach, see [10, 18]). 



the product of several factors, namely 

^l:^ = ^-tlc..e(.)x(^-)/^r--™(f,$) 

X ^i^^T,„.e'^-pL- (1) 

Here 6 = for "even-parity" (mass-generated) multi- 
poles {i + m even), and e = 1 for "odd-parity" (current- 
generated) ones {t + m odd); n^^, a-^d ci+e{i^) are numer- 
ical coefficients; 5*^^ is a /i-normalized effective source 
(whose definition comes from the EOB formalism) ; is 
a resummed version [1, 2] of an infinite number of "lead- 
ing logarithms" entering the tail effects [19, 20]; 5irn. is 
a supplementary phase (which corrects the phase effects 
not included in the complex tail factor Tim), and, finally, 
{.Pimf denotes the ^-th power of the quantity pira which 
is the new building block introduced and studied in this 
paper. [In previous papers [1, 2] the quantity {p(.rnf was 
denoted as fim-] 

We shall discuss in quantitative details below the var- 
ious facts showing that the new ingredient pim{x) is a 
useful quantity to consider (mainly because its PN expan- 
sion has better convergence properties than the straight- 
forward PN expansion of itself) . In this introductory 
section, we shall whet the appetite of the reader by com- 
paring the performance of our new resummed method, 
to some of the previously considered PN-based methods. 
For definiteness, we shall do this initial comparison at 
the level of the total energy flux, say F, which is related 
to the individual waveforms via 

t=2 m=l 1=1 m=l 

= i^EE("^")'i^^^'»i'- (2) 

1=2 m=l 

Note that Fim = ^£|m| denotes the sum of two equal 
contributions corresponding to +m and — m (m 7^ as 
FiQ vanishes for circular orbits). This explains the ex- 
plicit factor two in the last two equations above, which 
relate Fim to him- It is convenient to consider the total 
flux F for continuity with previous studies of the "con- 
vergence" of PN-expansions that focussed on F [4-6, 17] 
and because of its physical importance as a measure of 
the radiation reaction that acts on inspiralling binaries. 
To be fully precise, we shall consider here the (rather ac- 
curate) approximation i^^^-* obtained by truncating the 
sum over £ in Eq. (2) beyond £ = 6, and we normalize 
the result onto the "Newtonian" (i.e., quadrupolar) re- 
sult Fi^ = 32/5(/i/M)2x5. In other words, we consider 
here the quantity F = F^^) /F^. 

Fig. 1 compares and contrasts four different ways of 
using the same PN information about the total Newton- 
normalized GW flux function F{v) = F{v)/F22{v) (i.e., 
the same finite set of coefficients {fk{logx); 1 < k < n} 
of the n-PN expansion of the Newton-normalized flux 
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FIG. 1: Extreme-mass-ratio limit {1/ = 0). Comparing various resummations of the (Newton-normalized) gravitational wave 
energy flux: (a) standard Taylor expansion; (b) Pade resummation as proposed in Ref. [6] with lipoic = f/\/3; (c) Pade 
resummation flexing Vpoie according to the discussion of Sec. II of Ref. [2]; (d) new resummation technique based on the pirn 
functions discussed in this paper. 



Taylor„(/'(x)) = ELo /^(log^)^'' ^^e GW flux). As 
in many previous works, we use the extreme mass ratio 
limit 1/ ^ as a laboratory for devising and testing re- 
summation procedures. Indeed, in that case, the quasi- 
circular adiabatic description of inspiralling binaries is 
expected to hold up to the LSO (a;LSo('^ = 0) = 1/6) and 
one can compare PN-based analytical results [15, 16, 21] 
to numerical estimates of the GW flux, based on the 
numerical integration of some (Regge-Wheeler-Zerilli or 
Teukolsky) wave equation [4, 17]. 

Panel (a) of the figure recalls the results of Refs. [4, 5], 
namely the rather poor convergence of the standard Tay- 
lor approximants of F(x) in the full interval < x < xlso 
where one might hope to tap the PN information. For 
clarity, we selected only three Taylor approximants: 3PN 
(v^), 3.5PN (w^) and 5.5PN (z;"). These three values suf- 
fice (by contrast with the other panels) to illustrate the 



rather large scatter among Taylor approximants, and the 
fact that, near the LSO, the convergence towards the ex- 
act value (solid line) is rather slow, and non monotonic. 
[Sec Fig. 1 in Ref. [5] and Fig. 3 of Ref. [6] for fuller il- 
lustrations of the scattered and non monotonic way in 
which successive Taylor expansions approach the numer- 
ical result.] 

On the other hand, panel (b) recalls the results of [6], 
namely the significantly better (and monotonic) way in 
which successive Pade approximants approach (in L^c 
norm on the full interval < a; < a;Lso) the numerical 
result. Ref. [6] also showed that the observationally rele- 
vant overlaps (of both the "faithfulness" and the "effec- 
tualness" types) between analytical and numerical adia- 
batic signals were systematically better for Pade approx- 
imants than for Taylor ones. Note that our present panel 
(b) is slightly different from the corresponding results in 
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TABLE I: Errors in the flux of the two (untuned or tuned) 
Pade resummation procedures. From left to right, the 
columns report; the PN-order; the difference between the re- 
summed and the exact flux, AF - /'^-^"■"■"^d _ ^Exact^ 
the LSO, and the Loo norm of AF, ||AF||oo (computed over 
the interval < v < ulso), for Upoic ~ l/\/3; the flexed value 
of Wpoic used here; AF at the LSO and the corresponding Loo 
norm (over the same interval) for the flexed value of iipoic. 
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0.048 
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0.051 
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FIG. 2: Extreme-mass-ratio limit {v = 0). Complement to 
panel (d) of Fig. 1. Difference between the resummed and 
exact energy flux, for different approaches to the resummation 
of the P22 function. See text for explanations. 



panel (b) of Fig. 3 in [6] (in particular, the present panel 
(b) exhibits a better "convergence" of the v^-^ curve). 
This difference is due to the new treatment of the loga- 
rithmic terms oc log a;. Instead of factoring them out in 
front as proposed in [6], we consider them here (follow- 
ing [2]) as being part of the "Taylor coefficients" /n(log x) 
when Padcing the flux function. Note also that panel (b) 
follows Ref. [6] in introducing a pole in the resummed 

flux F{v) located at the value Wpo^°^ — l/VS- 

By contrast, panel (c) of the figure illustrates the re- 
markable improvement in the (Loo) closeness between 
^Padc-resummed^^) and F^xact^^) obtained, as recently 
suggested by Damour and Nagar [2] (following ideas orig- 
inally introduced in Ref. [22]), by suitably flexing the 
value of Vpoic- As proposed in Ref. [2], Vpoie is tuned 
until the difference between the resummed and the ex- 
act flux at the LSO is zero (or at least smaller than 



10~^). The resulting closeness between the exact and 
tuned-resummed fluxes is so good (compared to the pre- 
vious panels, where the differences were clearly visible) 
that we need to complement panel (c) of Fig. 1 with Ta- 
ble L This table compares in a quantitative way the re- 
sult of the "untuned" Fade resummation (upoic ~ ^/V^) 
of ReL [6] to the result of the "wpoic-tuned" Fade resum- 
mation of Ref. [2]. Defining the function AF{v;Vpoie) = 
^Rosummcd(-„. _ ^Exactj'^) measuring the difference 
between a resummed and the exact energy flux. Table I 
lists both the values of AF at v = wlso Sund its Foo norm 
on the interval < u < wlso for both the untuned and 
tuned cases. Note, in particular, how the Upoic-flcxing ap- 
proach permits to reduce the Laa norm over this interval 
by more than an order of magnitude with respect to the 
untuned case. Note that the closeness between the tuned 
flux and the exact one is remarkably good (4.3 x 10"'^) 
already at the 3PN level. 

Finally, panel (d) of Fig. 1 illustrates the even 
more remarkable improvement in the closeness between 
^Ncw-rosummod ^^^^ ^Exact obtained by means of the new 

resummation procedure proposed in this paper. More 
precisely, panel (d) plots two examples of fluxes obtained 
from our new p^m-representation, Eq. (1), for the individ- 
ual multipolar waveforms him in the sum Eq. (2) . These 
two examples differ in the choice of approximants for the 
£ = m = 2 partial wave. One example uses for P22 its 
3PN Taylor expansion, T3[p22], while the other one uses 
its 5PN Taylor expansion, r5[p22]- AH the other partial 
waves are given by their maximum known Taylor expan- 
sion. Note that the fact that we use here for the phn's 
some straightforward Taylor expansions does not mean 
that our new procedure is not a resummation technique. 
Indeed, the deflning resummation features of our proce- 
dure have four sources: (i) the factorization of the FN 
corrections to the waveforms into four different blocks, 
namely T<?„, e'"'*" and in Eq. (1); (ii) the fact 

the S^g is by itself a resummed source whose FN expan- 
sion would contain an infinite number of terms; (iii) the 
fact that the tail factor is a closed form expression (see 
Eq. (19) below) whose FN expansion also contains an 
infinite number of terms and (iv) the fact that we have 
replaced the Taylor expansion of fem = Pim by that of 
its i-th root, namely pim- 

Even more so than in the Wpoic-tuned case of panel 
(c), the closeness between analytical and exact results 
exhibited by the "new-resummed" case of panel (d) is so 
good that it is undistinguishable by eye. We therefore 
complement panel (d) by displaying in Fig. 2 the resid- 
ual differences AF{v) = /'Ncw-rcsummcdj-^) _ ^Exact^-^)^ 

We included in Fig. 2 a third curve corresponding to the 
case where we further resum our "new-resummed" fiux 
by using for (the 5FN accurate) P22 its near-diagonal 
(2,3) Fade approximant, say P| {T5[p22(a;)]}, instead 
of its Taylor expansion. [The other p^m's being still 
used in Taylor expanded form]. Note that the difference 
AF at the LSO is w 4.5 x 10"^ when using T3[p22], is 
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—5.7 X 10~^ when using T^[p22\ and is +1.6 x 10"'^ when 
using {T^[p22{x)W ■ Note that these numbers are in 
the same bah park than the f ^^-accurate Wpoio-tuned re- 
sult (8.5 X 10~^) quoted in Table I. Discarding the very 
small difference corresponding to the 3PN-accurate T-^ 
case as being probably accidental we conclude that us- 
ing the normal, near-diagonal^ Padc rcsummation of only 
the leading multipolar contribution p22 has the effect of 
significantly improving the agreement with the exact re- 
sult (compare the two 5PN-accurate curves, T^\p22\ and 
{T^[p22{x)]} ■, in Fig. 2). We therefore expect that 
Padeing some of the higher multipoles will further im- 
prove the agreement between the energy fluxes. Note 
also, in passing, that the new resummation procedure 
explored here is more "predictive" than the Upoic-flexing 
technique in that it does not need to rely on the knowl- 
edge of the exact answer. We will also show below that 
it is "robust" under the deformation brought about by 
increasing the symmetric mass ratio v from up to its 
maximal value 1/4. 

This paper is organized as follows: in Sec. II we ex- 
plicitly define the building blocks entering the resumma- 
tion of the gravitational waveform. Section III analyzes 
the performance of this resummation procedure in the 
extreme-mass-ratio [v 0) case, while Sec. IV consid- 
ers the comparable mass case (7^ 0). We present some 
concluding remarks in Sec. V. The paper is completed 
by three Appendices which collect many useful formulas 
that complete the discussion of the main text. In partic- 
ular. Appendix A provides the first explicit, ready-to- use 
expression for the 1-PN corrections to the Symmctric- 
Trace-Free (STF) current-multipolc moments and to the 
corresponding spherical-harmonics odd-parity multipoles 
for arbitrary i and m. 



five decomposition will be the Newtonian waveform, 
Some of the other factors are chosen so as to best capture 
part of the essential physics contained in the waveform. 
The remaining factors will then resum the subleading ef- 
fects that have not been included in the previous ones. 
First of all, it is convenient to introduce the following 
notation 

him = h'^fm"^ hf^, (3) 



where represents the Newtonian contribution and 

^hn product of all the other factors in our multi- 
plicative decomposition. As all these other factors repre- 
sent resummed version of PN corrections (of the type 
1 -I- 0{x)) their product will also have the structure 

h^'hm = l + 0{x). 



A. The Newtonian multipolar waveform 

Though all the work in this paper will focus on the 
resummation of the PN-correcting factor /i^^, let us, for 
completeness, recall the well-known [12, 14, 23] structure 
of the Newtonian multipolar waveform^, here considered 
for the adiabatic circular case. The Newtonian contribu- 
tion for circular orbits is, for given (£, to), a function of 
X = {GMVl/c')'^/^ and the two mass ratios Xi = mi/M 
and X2 = TOa/M^ 



II. DEFINING THE BUILDING BLOCKS OF 
THE RESUMMATION OF THE MULTIPOLAR 
GRAVITATIONAL WAVEFORM 

Let us now explicitly define each of the building blocks 
of our new resummed waveform, Eq. (1). Note that our 
methodology differs from the PN-methodology in a basic 
way. The PN approach consists in writing any relativistic 
quantity as a sum of various contributions starting with 
the so-called Newtonian approximation. In other words, 
a PN-expanded multipolar waveform has the structure 
hi-m = h^^ + h]™ + hl;f^ + .... By contrast to 
this additive approach we will use here, as advocated 
in Refs. [1, 2], a multiplicative approach in which any 
relativistic quantity is decomposed as a product of vari- 
ous contributions^. One of the factors of this multiplica- 



* We will explore other Pade approximants of p22 below. 
^ This multiplicative approach can be naturally applied to the mul- 
tipolar waveform which is a complex number. 



= ^42c...(^).^^-^/^l^— (f , . (4) 

Here, e denotes the parity of the multipolar waveform, 
i.e.; even (e — 0) for mass-generated multipoles and odd 
(e = 1) for current-generated ones. In the circular case, 
e is equal to the parity of the sum £ + m: e = tt{£ + to). 
In other words e = when ^ -|- to is even, and e = 1 when 
£+m is odd. The Y^^"-{9, 0) are the usual scalar spherical 



® We mostly follow here the conventions of Refs. [12, 23], except 
that we take into account some of the simplifications used in [14]. 
Note the presence of a factor 1 / %/2 in the relation between the 
{£, m) Newtonian waveform and the corresponding {£, m) radia- 
tive multipoles; = {G/V2)Uem = {G/V2)I^/2 ^'^^ 
multipoles and Rh°^ = -i(G/V2)Vf„ = -i{G/V2)5^2 for cur- 
rent multipoles. 

Note that: Xi + X2 = 1, X1X2 = v where v is the symmetric 
mass ratio v = m,\m2l{ra\ -|-r?i2)'^, while X\ — X2 = sign(mi — 
m2)\/l ~ 41^, where the sign depends whether mi > m.2 or the 
reverse. 
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harmonics while 



(0) ,. si> Stt 



(1) (■ \ 



(2£+l)!!y ^(£-1) 



(5) 



(6) 



are (^/-independent) numerical coefRcients. Fi- 
nally, the ;/-dcpendcnt coefficients Cf-|_c(;/) (such that 
I €£+£(;/ = 0)1 = 1), can be expressed in terms of v (as 
in Ref. [12, 14]), although it is more conveniently written 
in terms of the two mass ratios Xi and X2 in the form^ 



c,Uu)=Xf,+^-^ + {-Y+^Xi+^-' 



^2 



rn Y'^+€— 1 



(7) 



In the second form of the equation we have used the fact 
that, as e = 7r(£ -|- m), it{£ + e) = 7r(m). 



B. The first three factors in the multiplicative 
decomposition of the PN fractional correction /i^^ 

Let us recall that, in the comparable mass case {v 7^ 
0), the PN-correction can be computed within the 
perturbative multipolar-post-Minkowskian (MPM) for- 
malism [3, 19, 20], while in the test-mass limit {v — > 
0) it can be obtained by black hole perturbation the- 
ory [15, 16, 24, 25]. 

The final result is that lif^^ is given by a PN expansion 
of the form him = 1 + hix + hi^^x'^''^ + . . . . For 
comparable-mass circularized compact binaries, the par- 
tial wave which is known with the highest PN accuracy 
is the leading even-parity quadrupolar wave /122, which 
is known to fractional 3PN accuracy [2, 12, 14, 26]. 
Note that Ref. [14] provides him half a PN order more 
accurately than Ref. [12] for multipolar orders {£, to) = 
(2, 1), (3, 3), (3, 2), (3, 1), (4, 3), (4, 1), (5, 5), (5, 3) ^ and 
(5,1). This information is fully employed in this work. 
In the extreme-mass-ratio case, the partial waves are 
known with even higher PN accuracy. For instance, 
/i22 is known to 5.5PN [15, 16] and other multipoles to 



8 We use the Ye^'s defined in Eqs. (2.7) and (2.8) of Ref. [23], or 
equivalently by the s = case of Eqs. (4) and (5) of Ref. [12]. 

^ When expressing ci^i,{i') as an explicit function of u, as in 
Ref. [12], it is useful to note that cij^^{v) vanishes in the equal 
mass case when ^ -|- e is odd, which is equivalent (given that 
e = 7r(£ -I- m) for circular orbits), to m being odd. In such 
cases one can factor out, as in Ref. [12, 14], from cij^f^{u) a 
factor A = X\ — X2 = sign(mi — m,2)\/l — 4!/. The ci{u) 
in this paper are the same as the Si{v) in [14]. Note how- 
ever that in Appendix A of Ref. [12], in the second line above 
Eq. (A7), the definition of should include a supplementary 
factor m/&m = l/(Xi — X2) on the right hand side. 



accuracies consistent with 5.5PN GW flux. As explained 
later, this information is also appropriately exploited in 
our construction. 

As indicated above, the resummation method we shall 
use here consists in: (i) decomposing the PN-correction 

factor h!f^^ = 1 + hix + hi ^jX^^"^ + ■ ■ ■ into the product of 
four factors, each of which has a similar PN expansion, 
1 + 0{x), namely 



(8) 



and then (ii), resumming separately each factor. 

The choice of these various factors is based on our phys- 
ical intuition of the main physical effects entering the fi- 
nal waveform. The first factor is motivated by thinking 
about the form of the equation satisfied by each partial 
wave in the (circular) test-mass limit: indeed in this limit 
him is the asymptotic value (at spatial infinity) of a solu- 
tion of a (frequency-domain) wave equation of the Regge- 
Wheeler-Zerilli type (see e.g. Ref. [27]). The source term 
appearing on the right-hand side of this equation is a lin- 
ear combination of terms linear in the stress-energy ten- 
sor T^i, of a test-particle of mass fi moving around a black 
hole of mass M. As the Effective-One-Body method has 
shown that the dynamics of comparable-mass black holes 
can be mapped onto the dynamics of an effective particle 
of mass fi moving in some effective metric (which reduces 
to the Schwarzschild metric of mass M when j/ — > 0), it 
is natural to introduce (both when ly ^ and v ^ 0) 
effective source terms in the partial waves made up from 
the important dynamical characteristics of the FOB dy- 
namics, namely the effective FOB Hamiltonian Hcs and 
the FOB angular momentum J^. This motivates us to 
define as first factor in hfm an effective source term S, 



cff 



proportional either to Hctf or J'. Note that this idea of 
factoring Heft or J7 from the wave amplitude is similar 
to the suggestion of Ref. [6] of factoring out a pole in 
the energy flux. Indeed, the analytical continuation in x 
of the flux function F{x) below the LSO inherits, in the 
ly limit, a simple pole from the fact that F{x) is pro- 
portional to the square^° of the energy of the "rotating 
source" (see discussion p. 893 of [6]). 

Our second factor is motivated by thinking about 
the structure of the "transfer" function relating (in the 
comparable mass case) the far-zone GW amplitude him 
to the near zone one. If we keep, in the full Fin- 
stein equations considered outside the binary system, 
only the terms coupling the instantaneous "monopo- 
lar" Arnowitt-Deser-lVIisner (ADIvI) mass of the system, 
^^ADM = M + binding energy = //real, to the multi- 
polar wave amplitude, we get (in the circular approxi- 
mation and in the Fourier domain) a Schrodinger-type 



Note that in the u limit both -ffeff a^nd •J have a square-root 
singularity oc 1/Vl — 3x at the light-ring. Sec e.g. Eqs. (56)-(58) 
below. 
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equation, for each multipolc order containing a poten- 
tial Vi (r) whose leading behavior as r — *■ cxd is dominated 
by two effects: (i) the £{£ + l)/r^ centrifugal barrier, 
and (ii) a more slowly decreasing term ^ —AMabm^'^/t 
coming from the coupling to a curved (Schwarzschild- 
like) "monopolar" background metric. One can solve this 
leading-order equation by means of Coulomb wave func- 
tions. When doing this, it is found that each asymptotic 
partial wave is related to its corresponding near-zone ex- 
pression by a certain "tail" factor T^^. It can be checked 
that, in the comparable mass case, this tail factor rep- 
resents the resummation of the infinite number of lead- 
ing logarithms (see Eqs. (7)-(9) in [2]) that appear when 
computing asymptotic multipolar waves in the MPM for- 
malism [3, 19, 20]. Having so factorized two of the main 
physical effects entering /if™, we define the two other 
factors as the phase, e"^*'" and the modulus, fim of the 
remaining quoticnted Newton-normalized waveform. In 
this subsection we discuss in detail the first three fac- 
tors, postponing to the following subsection the last one, 
namely the modulus fim- 

Let us start by discussing the structure of the S'^g and 
Tim factors. In the even-parity case (corresponding to 
mass moments), since the leading order source of grav- 
itational radiation is given by the energy density, it is 
natural to define 



m even. 



(9) 



Here, ifoff is the effective EOB Hamiltonian (per unit 
fi mass), that we shall restrict here along the sequence 
of EOB circular orbits. When u 0, H^fi reduces 
to the usual conserved energy of a test-mass /i in a 
Schwarzschild background of mass M (see Eq. (56) be- 
low). 

The explicit expression of i?cff, along circular orbits, 
as a function of the frequency parameter x, cannot be 
written in closed form [7, 8]. However, it can be written 
in parametric form in terms of the EOB inverse radius 



parameter 
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1/r. More precisely, we have 



(circular orbits). 



(10) 

where u = 1/r, and where A{u){= — 5oo°'^*'™('')) 
crucial EOB radial potential and j ~ J /{jiGM) is the 
(dimensionlcss) angular momentum along circular orbits. 
Wc recall that the PN expansion of A{u) has the form 



^Taylo.(^) =1 - 2U + 2VU^ + (^y - 

+ 0(z^u^). 



(11) 



As usual in EOB work, wc use dimensionless variables, notably 
r = Rc^ /GM, where R is the EOB Schwarzschild-like radial 
coordinate. 



where the term corresponds to 3PN contributions to 
the EOB dynamics [28] and where we have parametrized 
the presence of yet uncalculated 4PN (and higher) contri- 
butions to A{u) by adding a term -\-as{v)u^ with the sim- 
ple form 05(1^) = a^v. As in previous EOB work, wc shall 
not use the Taylor-expanded function ^'^''^'"''(u), but re- 
place it by a suitably Pade resummed function A{u). 

The circular orbits in the EOB formalism are deter- 
mined by the condition du {A(m)[1 + pu^]} — 0, which 
leads to the following parametric representation of the 
squared angular momentum: 



(uM(u))' 



(circular orbits), (12) 



where the prime denotes d/du. Inserting this u- 
parametric representation of in Eq. (10) defines the 
M-parametric representation of the effective Hamiltonian 
Hcff{u). Wc can then obtain (at least numerically) Hcs 
as a function of x by eliminating u between H(.s{u) and 
the corresponding it-parametric representation of the fre- 
quency parameter x ~ {GM^l/ c?)'^/'^ obtained by the 
angular Hamilton equation of motion in the circular case 



Mn{u) 



1 dHi-cal 



MA{u)j{u)u'^ 



where i?rcai denotes the real EOB Hamiltonian 



real 



7\/Jl + 2z/ iJcff- 1 



(13) 



(14) 



While in the even-parity case we shall factor out 
Hcff{x) as a "source term", in the odd-parity one we ex- 
plored two, equally motivated, possibilities. The first one 
consists simply in still factoring Hcff{x); i.e., in defining 



= iJcff (x) £ + m odd. 



(15) 



The second one consists in factoring the angular mo- 
mentum J7. Indeed, the angular momentum density 
tijkX^T^^ enters as a factor in the (odd-parity) current 
moments, and occurs (in the small-i/ limit) as a fac- 
tor in the source of the Regge- Wheeler- Zerilli odd-parity 
multipolcs. This leads us to define as second possibility 



S'-^'^ ^K^)^x''^]{x) 



' + m odd. 



(16) 



where j denotes what can be called the "Newton- 
normalized" angular momentum, namely the ratio j(x) = 
j{x)/jN{x) with jNix) = \l\fx. [This Newtonian nor- 
malization being such that j{x) = 1-1- 0(x)\ We will 
compare below the performances of these two possible 
choices. Note that the PN expansions of these two pos- 
sible sources start as 



H.s{x) = 1 - + 0(a;2), 



(17) 



j-(.T) = l+(^ + f)x + 0(.T2). (18) 
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The second building block in our factorized decomposi- 
tion is the "tail factor" Tg^ (introduced in Refs. [1, 2]). 
As mentioned above, Tim is a resummed version of an in- 
finite numbers of "leading logarithms" entering the trans- 
fer function between the near-zone multipolar wave and 
the far-zone one, due to tail effects linked to its propa- 
gation in a Schwarzschild background of mass A/adm = 
H^^eid- Its explicit expression reads 



T{l+l- 2ik) 



Tk 2ik\og{2kro) 
^ 1 



(19) 



where = 2GM and k — Gi/rcai'Tifi and k = mQ. Note 

that k differs from fc by a rescaling involving the real 
(rather than the effective) EOB Hamiltonian, Eq. (14). 

The tail factor Tgm is a complex number which already 
takes into account some of the dcphasing of the partial 
waves as they propagate out from the near zone to infin- 
ity. However, as the tail factor only takes into account 
the leading logarithms, one needs to correct it by a com- 
plementary dephasing term, e"'*" , linked to subleading 
logarithms and other effects. This subleading phase cor- 
rection can be computed as being the phase Sgm of the 
complex ratio between the PN-expanded h'f^^ and the 
above defined source and tail factors. In the comparable- 
mass case {ly ^ 0), the 3PN S22 phase correction to the 
leading quadrupolar wave was computed in Ref. [2] (see 
also Ref. [1] for the ly ^ limit). For the subleading 
partial waves, we computed the other Sim's to the high- 
est possible PN-accuracy by starting from the currently 
known 3PN-accurate i/-dependent waveform [14]. Our 
explicit results read 



S22 = ly'^' 



S21 



"33 



2 3/2 _ 

3^ 

10 



4287r 

105 
493;^ 



S32 = 



10 + 33i/ 



15(1 - 3j/)" 



42 ^ 
80897 

2430 

3/2 



y - 24i/?/ 

5/2 



5/2 



5/2 



031 



043 



H„3/2 

30^ 
112 



To 
y 



y 

3/2 



5/2 



120(1 -3i/) 
486 + 4961:/ 
810(1 -2;^; 



-y 



3/2 



7(1 + 6z/) 3/2 
042 = -nrr^ — tttV : 



S41 = 



3/2 



055 



15(1 -3z/)" 
2 -I- 507v 
10(1 -2i/)^ ' 
96875 + 857528^ 
131250(1 - 2^) 



■y 



3/2 



(20) 
(21) 
(22) 
(23) 
(24) 
(25) 
(26) 
(27) 
(28) 
(29) 



Here, following Ref. [2], we define y = {H^cai^y^^ , which 
gathers together relativistic corrections (like those en- 
tering the tail) that depend on the instantaneous ADM 



mass of the system, namely i?roai, rather than the total 
"mechanical mass"M. Concerning the last y^^^ correc- 
tions^^, it is not clear whether they are more linked to 
the ADM mass or to the mechanical mass. This is why 
we use the notation y, meaning that it could be replaced 
either by x or y [note that in Ref. [2] we chose y = x 
inside the — 24^y^/^ correction to 522, Eq. (11) there ]. 
Indeed, these 2.5PN terms are not known to IPN frac- 
tional accuracy because we rely here on the available 3PN 
(and not 3.5PN) accurate results of [14]. 

In the extreme-mass-ratio limit ^ ^ 0, the informa- 
tion needed to compute some of the higher-order PN 
corrections to the iJim's is contained in the results of 
Ref. [15]. We leave to future work the task of exploit- 
ing this information to complete the above z/-dcpendcnt 
iSfm's with higher-order v = Q corrections. In addition 
we shall leave here the Sim's in Taylor-expanded form. 
We leave to future work an eventual comparison between 
numerically determined phases and (possibly resummed) 
analytic ones. 



C. The fourth factor in the multiplicative 
decomposition of the PN fractional correction /i 



(0 



The fourth and last factor in the multiplicative decom- 
position, Eq. (8), can be computed as being the modulus 

fim of the complex ratio between the PN-expanded /i^^ 
and the above defined source and tail factors. In the com- 
parable mass case (^ ^ 0), the /22 modulus correction to 
the leading quadrupolar wave was computed in Ref. [2] 
(see also Ref. [1] for the v = limit). For the subleading 
partial waves, we compute here the other fim's to the 
highest possible PN-accuracy by starting from the cur- 
rently known 3PN-accurate i/-dependent waveform [14]. 
In addition, as originally proposed in Ref. [2], to reach 
greater accuracy the fim{x; ^)'s extracted from the 3PN- 
accurate v results are complemented by adding 
higher order contributions coming from the j/ = re- 
sults [15, 16]. In the particular /22 case discussed in [2], 
this amounted to adding 4PN and 5PN 1/ = terms. This 
"hybridization" procedure is here systematically pursued 
for all the other multipoles, using the 5.5PN accurate 
calculation of the multipolar decomposition of the grav- 
itational wave energy flux done in Refs. [15, 16]. It is 
worth emphasizing at this stage that our hybridization 
procedure is not equivalent to the straightforward hybrid 
sum ansatz, him = h}m^''{v) + ^^f "('^ = 0) (where 
him = him/v) that one may have chosen to implement. 
The detailed definition of the hybridization procedure 



^'^ Note that these O (jj^^'^) corrections in arc the only terms in 
the wave ampUtude that we use here which go beyond the strict 
'circular limit' in that they include contributions proportional 
to the (radiation-reaction-driven) time derivatives of the orbital 
radius, or of the orbital frequency. 
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that we use, as well as the reasons why we think that 
our procedure is better than others, will be explained 
below. 

In the even-parity case, the determination of the mod- 
ulus fim is unique. In the odd-parity case, it depends 
on the choice of the source which, as explained above, 
can be either connected to the effective energy or to the 
angular momentum. We will consider both cases and dis- 
tinguish them by adding either the label H oi J oi the 
corresponding /£,„• Note, in passing, that, since in both 
cases the factorized effective source term (TJeff or J) is a 
real quantity, the phases S^rnS are the same. 

The above explained procedure defines the /^m's as 
Taylor-expanded PN series of the type 

u) = 1+4'^ {iy)x+4'^ {i,, \og{x))x'' + . . . 

(30) 

Note that one of the virtues of our factorization is to 



have separated the half-integer powers of x appearing in 
the usual PN-expansion of h^/^ from the integer powers, 
the tail factor, together with the complementary phase 
factor e"''^" , having absorbed all the half-integer powers. 

We have computed all the /fm's (both for the H and 
J' choices) up to the highest available (i^-dependent or 
not) PN accuracy. In the formulas for the fim^s given 
below we "hybridize" them by adding to the known v- 
dependent coefhcients c^'™ (i^) in Eq. (30) the i' ~ value 
of the higher order coefficients: cl^r'{v = 0). The IPN- 
accurate /£,„ 's for (. + m even and - thanks to the new 
results for IPN current multipoles obtained in Appendix 
A for arbitrary f - also for l+m odd can be written down 
for all i. In Appendix B we list the complete results for 
the ftmS that are known with an accuracy higher than 
IPN. Here, for illustrative purposes, we quote only the 
lowest ™" and fg^^"^ up to £ = 3 included. 



, , 1 , (2047t/2 _ 67451/ -4288) „ 
/22 (x; v) = l + -{hhv - 86)x + ^ — 

1146351/3 227875Z/2 41 . 34625i/ 856 , , , , 21428357\ . 

+ — TT 1/ — — — ■eulerlog2(.i;) H — 1 x 



99792 33264 96 3696 105 ' 727650 J 

, „o808 , , , , 5391582359\ , /458816 , , , , 93684531406\ 

culcrlognlx) a; + eulerlogn(x) 

V 2205 ' 198648450 J \ 19845 ^ 893918025 J 



, , , /23i/ 59\ /85i/2 269i/ 5\ , /88404893 214 , , , ^^ ^ 

fni (x;v) = 1+ \ ] x + [ a- + culcrlog, ix) x'^ 

■'^^^ ' ' V 42 28/ V 252 126 9/ V11642400 105 ') 

/6313 , , , , 33998136553\ 4 5, 

+ culerlog, (x) x^ + O(x^), (32) 

Vl470 ' 4237833600 / ^ ^ 1^ v 7 

^ , , , / 7\ /8871/2 MQlv 443\ 2 /147471561 78 , , , A 

t?,?,\x\v) = 1 + Zv .T + a; + eulerlogQ x x: 

■'^^^ ' ' ^ 2/ V 330 330 440y V 2802800 7 ^^"^ n 



39 eulerlog3(x) - ^f^) + 0{x% (33) 



,T , , , 3201/2 - 1115;/ + 328 39544i/3 _ 253768z/2 117215i/ - 20496 , 

fo,(x: I/) = 1 H ^ X H ; ^ X 

■'^^^ ' ' 90(31/ - 1) 11880(3;/ -1) 

110842222 104 , , 



4729725 21 

^ , , , / 2^ 13\ / 247!/- o... ..,0 , 2 

/400427563 26 . , , A , /169 , , , , 12064573043\ 



1 j x^ H 


-0{x^] 


3711/ 


1273 


198 


^ 792 


]x^ + 


/169 



V 75675600 21 ' ) V 63 ' 1816214400 / 



(34) 



x-^ + 0{x^). (35) 



For convenience and readability, we have introduced the following "eulerlog" functions eulerlog„j(a;) 

eulerlog„(a;) = 7^ + log 2 + ^ log .x + log to, (36) 
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which expHcitly reads, when m = 1, 2, 3, 



eulerlogi(x) = 7^ + log 2 + - logx, 
eulerlog2 (.t) = 7^ + 2 log 2 + i log 



(37) 
(38) 



eulerlogg (.t) = 7^ + log 2 + log 3 + - log x, (39) 

where 7e is Euler's constant, 7e = 0.577215 . . . and 
log(a;) denotes, as everywhere else in this paper, the nat- 
ural logarithm function. 



D. Resumming the modulus factor fim 

The decomposition of the total PN-corrcction factor 
several factors is in itself a resummation pro- 
cedure which has already improved the convergence of 
the PN series one has to deal with: indeed, one can see 
that the coefficients entering increasing powers of x in 
the fgm^ tend to be systematically smaller than the co- 
efficients appearing in the usual PN expansion of h!f^. 
The reason for this is essentially twofold: (i) the factor- 
ization of Tim has absorbed powers of mir which con- 
tributed to making large coefficients in h!f^, and (ii) the 
factorization of either H^s or j has (in the v ^ Q case) 
removed the presence of an inverse square-root singu- 
larity located at x = 1/3 which caused the coefficient 
of x" in any PN-expandcd quantity to grow as 3" as 
n — > cx). To prevent some potential misunderstandings, 
let us emphasize that we are talking here about a singu- 
larity entering the analytic continuation (to larger values 
of x) of a mathematical function h{x) defined (for small 
values of x) by considering the formal adiabatic circular 
limit. The point is that, in the v ^ Q limit, the ra- 
dius of convergence and therefore the growth with n of 
the PN coefficients of h{x) (Taylor-expanded at x = 0), 
are linked to the singularity of the analytically contin- 
ued h{x) which is nearest to x = in the complex x- 
plane. In the v case, the nearest singularity in the 
complex x-plane comes from the source factor Hcfi(x) 
or j(x) in the waveform and is located at the light-ring 
xuilv = 0) = 1/3. In the 1/ 7^ case, the EOB for- 
malism transforms the latter (inverse square-root) singu- 
larity into a more complicated ( "branching" ) singularity 
where dH^s/dx and dj/dx have inverse square-root sin- 
gularities located at what is called [2, 8, 11, 18, 29] the 
(Effective) ^"^ "EOB-light-ring" , i.e., the (adiabatic) max- 
imum of n, x%f^ii^) = {Mn^l^^f^^ > 1/3. 



Beware that this "EfTcctivc EOB-hght-ring" occurs for a circular- 
orbit radius slightly larger than the purely dynamical (circular) 
EOB-light-ring (where H^g and J' would formally become infi- 
nite) . 



Despite this improvement, the resulting "convergence" 
of the usual Taylor-expanded fimixYs quoted above does 
not seem to be good enough, especially near or below the 
LSO, in view of the high-accuracy needed to define grav- 
itational wave templates. For this reason, Refs. [1, 2] 
proposed to further resum the /22 (x) function via a Pade 
(3,2) approximant, P2 {/22(x; z^)}, so as to improve its 
behavior in the strong-field-fast-motion regime. Such a 
resummation gave an excellent agreement with numeri- 
cally computed waveforms, near the end of the inspiral 
and during the beginning of the plunge, for different mass 
ratios [1, 10, 18]. Here, however, we wish to explore a new 
route for resumming ftm, based on replacing ftm by its 
^-th root, say 



(40) 



Our basic motivation for replacing fim by pim is the 
following: the leading "Newtonian-level" contribution 
to the waveform /i^^^^ contains a factor cu^rf^^^^-^v'^ where 
7'harm IS the hamionic radial coordinate used in the MPM 
formalism [30, 31] . When computing the PN expan- 
sion of this factor one has to insert the PN expansion 
of the (dimcnsionlcss) harmonic radial coordinate rharm, 
Tharm = x~^(l 4- Cix -f- 0(x^)), as a fuuctiou of the 
gauge-independent frequency parameter x. The PN re- 
expansion of [riiarm(2^)]^ then generates terms of the type 
x^^{l + £cix + ....). This is one (though not the only one) 
of the origins of IPN corrections in and fim whose 
coefficients grow linearly with £. As we shall see in de- 
tail below, these i!-growing terms are problematic for the 
accuracy of the PN-cxpansions. Our replacement of fim 
by Pirn is a cure for this problem. 

More explicitly, the investigation of IPN corrections to 
GW amplitudes [12, 30, 31] has shown that, in the even- 
parity case (but see also Appendix A for the odd-parity 
case). 



_ Ci^2{v) m2(£ + 9) 
Ci(y) 2(^-hl)(2^-t-3)' 



1 ^ 3 CiJ^2{v) _ biji^) 

Ci{v) Ci{v) 



(41) 



where Ci{v) is defined in Eq. (7) and, consistently with 
the notation of Appendix A, 



bii,.) EE + {-Yx( 



(42) 



As we shall see below, the v dependence of c{''"' (v) is 
quite mild. For simplicity, let us focus on the ;/ = case, 
where the above result shows that the PN expansion of 
fim starts as 

(43) 

The crucial thing to note in this result is that as £ gets 
large (keeping in mind that |m| < £), the coefficient of 
X will be negative and will approximately range between 
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-5£/4 and -£. This means that when £ > 6 the IPN 
correction in /f„i would by itself make fim{x) vanish be- 
fore the {v = 0) LSO X — 1/6. For example, for the 
^ = ?Ti = 6 mode, one has /g™(a;;0) = 1 — 6a;(l + 
11/42) « 1 — 6.t(1 + 0.26) which means a correction equal 
to -100% at a; = 1/7.57 and larger than -100% at the 
LSO, namely /6^6^^(l/6; 0) w 1 - 1.26 = -0.26. This is 
definitely incompatible with the numerical data we shall 
quote below. Similar results hold also for the odd-parity 
fim^s, especially in the case where we factorize the J 
source (which happens to have close similarities to the 
iJ- factored even-parity fim)- Indeed, we have extended 
the result of Eq. (43) to the odd-parity case, i.e. we have 
computed, (using the comparable mass IPN results of 
Ref. [31]) the IPN correction in fgm and pe-m- In the 
v limit, we found that (see Appendix A for more 
details and for a discussion of the comparable-mass case) 



2 m^ii + A) 



(44) 



which is structurally similar to the even-parity expression 
quoted above. 

Let us now see how the replacement of fim by the 
newly defined pim, Eq- (40), cures this problem of abnor- 
mally large IPN corrections to the waveforms for large 
values of £. Indeed, the Taylor expansion of pim now 
starts as (say for simplicity in the v = 0, even-parity 
case) 



where, for £ ^ I, the coefficient of x again approximately 
ranges between —5/4 and —1. 

We have computed all the pim^s (both for the H and 
J choices) up to the highest available (z/-dependent or 
not) PN accuracy. In the formulas for the p^m's given 
below we "hybridize" them by adding to the known ly- 
dcpcndent coefficients cj^*™ {i>) in the Taylor expansion of 
Pern's, 



pira{x\v) = H-cf'"(i/)x-f-4'''"(i/)a;+cg*"(log(x);i^)a;^-|-... 

(47) 

the V = Q value of the higher order coefficients 
c''^r{v = Q)- Beware that this definition of an hybrid 
pgm is not equivalent to that displayed in Eqs (31)-(35) 
above of an analogous hybrid fpm (nor is it equivalent to 
a straightforward hybridization of him)- The primary hy- 
bridization procedure that we advocate (and use) in this 
paper is the one based on pim (i.e., replacing c^^'"(i') by 
c^'™ (0) when ??' is beyond the maximal j/-dependent PN 
knowledge). The IPN-accurate pim's for ^-f to even and 
- thanks to the new results for him for £ + m odd in Ap- 
pendix A - also for £ + m odd are explicitly known for all 
£. For the IPN coefficient of the p^m's we explicitly have 



_ z/\ 1 3 Ci+2{v) _ 1 hi{v) 
3/ 2£ 2£ Ci{v) £ci{v) 
m2(^ + 9) Ci+2{v) 



2£(£ + l)(2^-H3) Ci{v) 



(48) 



p,— (x;0) = l-x 1-- + 



TO^(^ + 9) 

2£{i + l){2£ + i) 



(45) 

Note that for large £ and arbitrary to the coefficient of 
X now approximately ranges between —5/4 and —1. We 
shall see below that the nice behavior of pim expected 
from this IPN estimate indeed holds for the exact pim., 
at least in the 1^ = case. In addition, the same struc- 
ture is found in the odd-parity p/^'s. In particular, from 
Eq. (44) above one finds 

p/m(2;;0) ^1-x 



1 _2_ m^{e + A) 

£ " £2 + 2e{£ + 2){2£ + 3) 



+ 0(x2), 
(46) 



2£- 



--)-- 
3/ 2£ 

3be+i{iy) , 



2£2 ci+i{iy) 

1£+1q+3(^) 
2 £2 c,+i(j.) 



2iy 



" 3/ 2£2 
£+lbf_i(t/) 

£2 Q+l(l.) 

m^{£ + 'i) 



2£{£ + 2){2£ + 3) Ci+i{v) 
(49) 



For definiteness, we give in Appendix B the complete 
results, for p^m (even-parity) and p/„^ (odd-parity), up to 
£ = 8 included. Here, for illustrative purposes, we quote 
only some of the lowest multipolc results up to £ = 3 
included. 
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, , :hhv 43\ /19583i^2 330251/ 20555\ , 



84 42/ V 42336 21168 10584 ) 

/10620745Z/3 6292061J/2 M-k'^v 48993925;/ 428 1556919113\ 3 

V 39118464 3259872 ^ ^^2 9779616 ^ 105'^'^ + 122245200 j 

/9202 , , , ^ 387216563023\ 4 /439877 , , , , 16094530514677\ ^ 



V2205 ' 160190110080/ V 55566 ' 533967033600 / 



(50) 



piiix;!') = 1 + 



23u 59 
"sT " 56 



a:; - 



/ 617^2 10993^ 47009 \ 
V 4704 " 14112 ~ 56448 j 



,2 



/7613184941 107 , , , x\ 3 /6313 , , , , 1168617463883\ 4 
+ (, 2607897600 " Wb'""''''"^^^^^ ) ^ + i5880^^^^^^°Si(x) - ,,,,,,,,,,,, ) + 0{x% (51) 

, , , {2v 7\ /149^2 igQi^ Qjig\ /3203101567 26 , , . ^^ 1 
P33(x; .) = 1 + (^y - - j X + j + 227026800 " Y^^l^'-l^&l-) j - 

/^13 , , , , 57566572157\ 4 5, 
+ -culcrlog3(.) - 33^2153600 ^ + )' ^^^^ 



J , , , 320i/2 - 1115J/ + 328 30856401/-* - 20338960i/3 - 4725605i/2 + 8050045;/ - 1444528 

P^nfx; V) = I -\ : ^ X H ; r;^ X 

^^^^ ' ' 270(3^- 1) 1603800(1- 3^)2 



2 



/ 5849948554 104 4, 



, , , , _Lo , , 829;/2 1685j/ 101 \ , / 11706720301 26 , , . ^^ 3 

P3i(x: i^) = 1 + he + \ a- + culcrlogi (x) a; 

'^'^^ ' ^ " 10) ;^782 1782 7128/ V 6129723600 63 7 



169 , , . , 2606097992581\ 4 5, 

culcrlog, (x) H x^ + 0(3;^). (54) 

567 '^ 4854741091200/ ^ > ^ ^ 



III. RESULTS FOR THE 
EXTREME-MASS-RATIO CASE {v = 0) 

A. Extracting the pf^'^^^ multipoles from black-hole 
perturbation numerical data 

To test our new resummation procedure based on the 
pim^s we shall compare the analytical results defined by 
our multiplicative decomposition Eq. (8) to the "exact" 
results obtained by numerical analysis of black hole per- 
turbation theory. For most of the comparisons discussed 
below we will rely on data kindly provided by Emanuele 
Berti, who computed the multipolar decomposition of the 
GW flux from stable circular orbits above the LSO with 
a frequency-domain code which solves numerically the 
Teukolsky equation with a point-particle source, (see for 
example Ref. [17] and references therein). In addition, we 



have complemented Fig. 5 by computing the quadrupolar 
GW energy flux F22 from a sample of unstable circular 
orbits with radius between 6Af and 3.1M. 

The result of the numerical computation is expressed in 
terms of the multipolar pieces Fim of the total "exact" 
flux, Eq. (2). We shall only consider multipoles up to 
£ ~ 6 included. The "exact" 1/ ^ version of our new 
quantities p^m(a;;^')'s are then obtained from the ratio 
of the "exact" partial fluxes F^^'^^ to their Newtonian 
counterparts -F^^**"" as 
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FIG. 3: Extreme-mass-ratio limit {v = 0). The "exact" functions pem{x) for < a: < 1/6 extracted from E. Berti's numerical 
fluxes. Multipoles up to ^ = 6 are considered. Each panel corresponds to one value of i and shows the even-parity partial 
amplitudes (black online) together with the J -normalized odd-parity ones (red online). 



where, for v = we explicitly have and either 
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VI - 3a; 



Still in the v ^ limit, we have ifrcai — >■ M, and, from 
well-known properties of the F function, the square mod- 
ulus of the tail factor T^m reads 



\Ten 



1 



47rfc 



1 - 



n + (2^-)' 



(59) 



where now k = mAIQ = mx^^"^ = 



B. Finding structure in the p^'''^* multipoles 
extracted from numerical data 



Let us first consider the properties of the "exact" pem's. 
In the odd-parity case we shall focus here on the J7- 
normalized quantity We shall see that the pernS 

convey interesting information about the x dependence 
of the multipolar GW amplitudes. Fig. 3 exhibits the 
numerical pim functions for i < 6, versus x = (MH,)^^^ 
(where we recall that fl denotes the orbital frequency) 
up to the I' = LSO, x ~ 1/6. Each panel of the figure 
displays, for each given 2 < £ < 6, the partial ptm's for 
the various possible m's, 1 < m < £ (we do not plot the 
negative m's since they correspond to the same value of 
Pirn)- The even-parity (£+m even, black online) and odd- 
parity {£ + m odd, red online) modes are shown together 
for comparison. 

Fig. 3 displays the following noticeable facts: (i) to a 
good approximation, all the ptm{x) are straight lines^'^ 
(see below) ; (ii) for each value of £, the (negative) slopes 
of the dominant m = £ (even-parity), and subdominant 
m ~ £—1 (odd-parity) multipole modes are very close to 
each other; and these slopes become closer and closer as 
the value of £ increases (note in particular that ior £ — 6 
Peei^) and pes (2;) are practically coincident); (iii) for a 
given value of £, and a given parity (even or odd), the 
absolute value of the (negative) slope decreases mono- 
tonically as |m| decreases. This "order" in the "exact" 
data can be analytically understood. 

The property (i) means that the IPN correction is al- 
ready capturing most of the physical information, which 
might turn out to be a useful fact to know (see below). 
We illustrate this result in Fig 4 which focuses on the 
quadrupolar {£ = 2) partial waves, and exhibits the ex- 
act P22 and P21 (solid and dashed lines) together with 



In the odd-parity case, i+m odd, this quasi-linear behavior up to 
the LSO, is particularly clear for the functions p^^(a:)'s, making 
us consider as our best default choice. By contrast, the 

i?-normalized functions p^^{x) have a more ^-dependent shape 
that the reader can figure out by noting the link between them: 



their IPN approximations (dotted lines). Note, for in- 
stance, that the difference between the IPN accurate. 



(58) Pa™, that we shall denote^^ Ti[p22] = P2™ = 1 



X, 

and the "exact" one Pa™'^* is equal, at the LSO, to 
P2™ - P2r''^ = 0.8294 - 0.8143 = 0.0151, which is only 
1.9% of the exact result 0.8143. The other multipoles 
exhibit a similar agreement between the exact pgm and 
their analytical IPN representations. To understand an- 
alytically what underlies this agreement, let us consider 
the £ = m = 2 case. Numerically, from Eq. (50) we have, 
near the LSO (for simplicity, we replace the log(a;) terms 
present in the coefficients by their numerical values at 
X = 1/6) 

P22{x] 0) « 1 - 1.024j; - lM2x^ 

+ 8.384a;3 -f 2.038a;^ - 21.690a;'^ 

« 1 - 0.171(6x) - 0.054(6a;)2 

+ 0.039(6a;)3 + 0.0016(6x)'' - 0.0028(6x)^ 

(60) 

We see that the successive coefficients of the PN- 
expansion of P22 are such that, even at the LSO, the 
magnitudes of the PN-corrections beyond the IPN one 
are rather small. They are significantly smaller than 
the corresponding terms in the usual PN-expansion of 
the total flux. For instance, by contrast to the coeffi- 
cient —21.69 which enters the 5PN-correction in P22, let 
us recall that the coefficient (including the log(a;) esti- 
mated at the LSO) of the 5PN correction in the usual 
PN-expanded flux is « -1321.402 (see e.g. [6]). Note 
that the latter 5PN contribution to the PN-expanded flux 
considered at the LSO is -1321.40/6^ « -0.17 which is 
as large as the IPN contribution to P22 and about sixty 
times larger than the corresponding 5PN correction to 
P22. In addition, as the signs in Eq. (60) fluctuate, there 
are compensations between the higher PN contributions, 
as it will be clear from further results presented below. 

Property (ii) can be analytically understood by means 
of the IPN-accurate closed formulas, Eqs. (45)-(46). In- 
deed, it is easily checked that the difference between the 
coefficients of x Eq. (45) for m = £ and Eq. (46) for 
TO = € — 1 is of order 0{l/£'^) when £ gets large. 

Finally, property (iii) is understood by noting that 
the coefficients of x in Eqs. (45)-(46) have the structure 
— {a{£) + m'^h{£)) where a{£) and h{£) are positive. 



C. Comparing Taylor and Pade approximants of p22 

Let us now compare and contrast the "convergence" of 
various PN-approximants towards the "exact" (numer- 
ical) P22. We first focus on the values of various ap- 



Ptm 



:pi„/(l-2x)l/^. 



Here and in the following we shall denote the truncated n-PN- 
accurate expansion of any function f(x) as Tn[f(x)] = /o-l-/l2;-l- 
■■■-h/„x". 
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FIG. 4: Extreme-mass-ratio limit {1/ — 0). Comparison be- 
tween the "exact" leading and subleading quadrupolar ampli- 
tudes P22 and P21 and the corresponding IPN-accurate ana- 
lytical ones. 



TABLE II: Closeness of various resummed approximants to 
p22 at the LSO, 2;lso = 1/6, or actually iiaat = 1/6.00001. 
The rightmost column lists the difference Ap22 between the 
resummed approximant and the exact value at a; = a;iast 



Approximant 


P22(2;iast) 


Ap22(a;iast) 


„Exact 

P22 


0.8143372247 





Ti[p22] 


0.8293653638 


0.0150281391 


T2[p22] 


0.7754188106 


-0.0389184141 


r3[p22] 


0.8142342355 


-0.0001029892 


^4^22] 


0.8158069452 


0.0014697205 


r5[p22] 


0.8130176477 


-0.0013195770 


Pt{TB[p22]} 


0.8148012716 


0.0004640469 


PHT5[P22]} 


0.8146954164 


0.0003581917 


Pi{TB[p22]} 


0.8132320684 


-0.0011051563 


PimP22]} 


0.8146954104 


0.0003581857 



proxinmnts to P22{x) at the LSO, i.e. at xlso = 1/6 or 
actually the last point in the numerical data computed 
by E. Berti, xust = 1/6.00001. At the point x = xust 
the numerical value of the Newton-normalized £ = m = 2 
partial flux is F22 = ^22/^^2^ = 0.8927266028. This cor- 
responds to pf^'^'^ixii.stj == 0.8143372247. In Table II wc 
compare this value to several PN-based approximants: 
both Taylor approximants, from IPN to 5PN (Ti[p22] 
to T5[p22]) and several of the "around the diagonal" 
5PN-accurate Pade approximants, namely, Pi{Tr,[p22]}, 
PI{T5[P22]}, Pi{T5[p22]} and P|{T5[p22]}. Note how the 
sequence of Taylor approximants to P22 nicely approaches 
the exact value, especially starting with the 3PN approx- 
imation. Probably by accident, the Taylor 3PN approx- 



imant. r3[p22]; happens to be closer to the exact value 
than the higher order approximants r4 and T^. Besides 
this accidental closeness of T3, the important thing to 
note is the very small dispersion (within ±1.8 x 10"'^) 
of Ts, Ti and around the correct value. This excel- 
lent behavior of the Taylor approximants of P22 should be 
contrasted with the much worse behavior of the standard 
Taylor approximants cither of the flux or of the wave- 
forms (see for example Fig. 1 in Ref. [5], Fig. 3 in Ref. [6] 
and Fig. 6 below, where we directly compare the usual 
Taylor approximants to the waveform to our new pem- 
based approximants). Note that when considering "Tay- 
lor approximants to P22" we are actually speaking of a 
specifically resummed approximant to the waveform /i22- 
This approximant has the factorizcd form of Eq. (8), and 
is made of the product of several resummed constituents. 
Even the last factor /22 of this product is not used in 
Taylor-expanded form (which would be T5[/22]), but in 
the minimally resummed way jRosummcd _ (X'5[p22])^. 

We have also explored several ways of further resum- 
ming P22, i-e., of replacing its PN-expanded form r5[p22] 
by various non- Taylor approximants. In view of the good 
closeness of the IPN approximation to P22 to the exact re- 
sult we explored, in particular, some "factorized" approx- 
imants (similar to those considered for the A{u) func- 
tion in Ref. [32]) of the type p22(a;) = (1 + cf "a;)p22(a;). 
We will not show our results for these approximants 
here. Instead, let us discuss the use of Pade approxi- 
mants for representing T^[p22\ as a rational function^® 
of X. As an example, we present in Table II the val- 
ues of p22(a;iast) predicted by using the four "around 
the diagonal" 5PN accurate Pade approximants, namely 

Pt{n[p22]}, Pl{n[p22]}, PiiTM) and Pl{n[p22]Y 
The important thing to note is that all these approxi- 
mants are both consistently clustered among themselves, 
as well as closely centered around the correct numerical 
value (within ±1.3x 10^"^). Note also that, apart from Ts, 
all the Pade approximants are closer to the exact value 
than the T/^ and T5. [Though the a priori less-accurate 
T3 approximant happens to be closer to the exact value 
than all other approximants, we consider that this is co- 
incidental because the subsequent Taylor approximants 
T4 and Ts do not exhibit such a close proximity.] 

Wc display in Fig. 5, the various 5PN-accurate approx- 
imants discussed above (together with the IPN-accurate 
Ti[p22] for comparison) to P22{x) over the larger inter- 
val Q < X < xlr, where xlr = 1/3 is the value of the 
frequency parameter x at the {v = 0) "light-ring" . Note 
that, while all the 5PN approximants stay very close to 
each other (and to the "exact" numerical value, red on- 
line) up to the LSO, they start diverging from each other 
when X > 0.2. This motivated us to extend the nu- 



As in Ref. [2] wo consider in this work the logarithmic terms (of 
the type eulerlog^ (a;) in Eqs. (50)-(54)) as part of the coefficients 
when Padeing pim{^) 
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TABLE III: Newton-normalized energy flux, and partial am- 
plitudes /22 and P22, for a sample of unstable circular orbits 
computed via the time-domain code of Ref. [33]. These val- 
ues of P22 are represented as empty circles in Fig. 5. The 
case r = 6 is shown here only for comparison with frequency- 
domain-based results. 



r 


X 


TTitimc 
-^22 


/■time 
722 


„timc 

P22 


6 


0.1666 


0.897 


0.665 


0.815 


5 


0.2000 


0.995 


0.615 


0.784 


4 


0.2500 


1.378 


0.562 


0.750 


3.5 


0.2857 


2.202 


0.539 


0.734 


3.1 


0.3226 


6.665 


0.513 


0.716 



mcrical data of E. Bcrti beyond the LSO, i.e. to con- 
sider the GW flux emitted by unstable circular orbits of 
Schwarzschild radii 3GM < R = GMr < 6GM, corre- 
sponding to 1/6 < x < 1/3. See Table III for results 
obtained for such a sample of sub-LSO orbits (they also 
appear in Fig. 5 as empty circles). These numbers have 
been computed with the time-domain code described in 
Ref. [33]. A resolution of Ar» = 0.01 was used. To test 
the accuracy of our numerical procedure we computed the 
energy flux F22 at r = 6 and compared it with the value 
obtained via Berti's frequency domain code (at the very 
close value r = 6.00001). We obtained ^2*^"° = 0.897342 
to be contrasted with i^22°'^ = 0.892726, which yields a 
fractional difference AF/F^'^'"i w 0.005. This gives an 
indication of the accuracy of our time-domain results, 
though we expect, for various numerical reasons, that 
the accuracy degrades as r gets below 4. 

We do not wish to give too much weight to the indi- 
cation given by our sub-LSO results on the behavior of 
the function P22{x) below the LSO. Indeed, on the one 
hand, the GW flux along sub-LSO circular orbits does 
represent the analytic continuation of the function F{x) 
yielding the GW flux along (stable) super-LSO circular 
orbits. As such, the empty circles in Fig. 5 do provide 
correct mathematical information about the analytical 
continuation of the function P22 (x) that we are trying to 
best approximate. On the other hand, we are evidently 
aware that the real dynamics of the "plunge" strongly 
deviates from the sequence of unstable circular orbits be- 
low the LSO and that the GW flux emitted by a plung- 
ing test-mass (or effective source) will not be correctly 
represented by this mathematical continuation of F{x). 
However, we expect, especially on the basis of the EOB 
formalism which has shown that plunging orbits remain 
approximately quasi-circular, that, in view of the present 
approach where we decompose the GW amplitude into 
several different factors having different physical origin, 
the mathematical continuation of the P22 part is likely 
to continue to capture important aspects of the nonlin- 
ear relativistic corrections to the waveform [note that we 
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FIG. 5: Extreme-mass-ratio limit (i^ = 0). Resummation of 
the function p22{x) on the interval < a; < 1/3: compari- 
son between various Taylor and Fade approximants and the 
"exact" function obtained from (both frequency-domain and 
time-domain) numerical calculations. The time-domain data 
points (see Table III) are indicated as empty circles. 



have in mind using our factorized waveform Eq. (8) along 
the EOB quasi-circular plunge together with the correct 
instantaneous source S'f^ and tail corrections Ti„i\ ■ How- 
ever, we are also aware that some aspects of the EOB 
plunge do physically differ, near the end of the plunge, 
in a relevant way from the physics included in the math- 
ematical continuation of p22{x): namely, the fact that 
the ratio {mViy' /Vi{r) (where Vt{r) is the Zerilli poten- 
tial) stays always small along the real plunge, while it 
increases more along unstable circular orbits and ends 
up reaching values of order unity. In other words the 
part of P22{x) which takes into account the filtering of 
Vf(r) will be different in the two cases for orbits near the 
light-ring. However, with due reserve we think that the 
first three empty circles on Fig. 5 do provide a guide- 
line for selecting among the various diverging PN ap- 
proximants the ones which are likely to provide, within 
the EOB formalism, a good zeroth-order approximation 
to the wave-amplitude emitted by real plunging orbits. 
But, we expect that it will be necessary to correct such a 
zeroth-order quasi-circular wave-amplitude by non-quasi- 
circular corrections of the type which has already been 
found necessary in Refs. [1, 18] to obtain a close agree- 
ment between EOB waveforms and numerical waveforms. 

If we use such a guideline, Fig. 5 suggests that the 
best continuations of ^22(2^) below the LSO are given by 
the three particular Fade approximants, Pi, and P|. 
However, as Pj* develops a spurious pole (which is barely 
visible on the left upper corner of the figure because it 
is very localized) at x w 0.038 we will discard it. By 
contrast, the other two are robust against the presence 
of spurious poles in the useful regime x < 1/3 (although 
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FIG. 6: Extreme-mass-ratio limit {y — 0): Various representations of the |/i22| waveform modulus. Left panels: standard 
PN-expanded amplitudes. Right panels: Various resummed amplitudes. See text for details. 



they develop poles for higher values of x, namely below 
the formal "event horizon" value x = 1/2). In the follow- 
ing, we shall choose P| as our current best-bet approx- 
imant to the p22 function (notably because this is the 
natural near-diagonal default Pade approximant). Note 
finally, in Fig. 5, how the simple IPN-accurate Taylor ap- 
proximant of P22{x) succeeds in providing a reasonably 
good representation of /9f^'^'^*(a;) over a very large range 
of X values. 



D. Comparing resummed waveforms to 
Taylor-expanded and exact ones 

Up to now, we focussed on the "convergence" of various 
PN-based approximants towards the numerically deter- 
mined value of the fourth technical building block P22 
entering the dominant quadrupolar wave. 

In this subsection we shall investigate instead the "con- 
vergence" of various possible PN-based approximants to- 



wards the more physically relevant Newton-normalized 
GW amplitudes hem- On the one hand, we shall con- 
sider not only the dominant i = m = 2 wave, but 
also a selection of subdominant partial waves. On the 
other hand, we shall consider other PN-based approxi- 
mants than those considered above. In particular, we 
shall compare and contrast the exact moduli \him\ both 
with standard high-accuracy Taylor-expanded waveforms 
(l^fml = 1 + cix + ci.scc^/^ -I- C2X^ -|- . . . ), and with our 
new "resummed with Taylor [p]" waveforms ( |/i£,„| = 

Scs\Tem\Pem ^^^^ Pim = 1 + c[x + c'^x'^ + Cga;^ -I- . . . ). 
We shall also analyze the performance of our new "re- 
summed with Pade[p]" waveforms ( \hi„i\ = S'off |T^m|p^m 
with ptm = Pp[l + c'lX + c'2X^ -t- CgX^ -t- . . .]), at least 
for the £ = m = 2 dominant mode.^^ For definiteness. 



In view of the remarkable agreement, displayed in panel (d) of 
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FIG. 7: Extreme-mass-ratio limit {u — 0). Resummation 
of the 1^221 waveform modulus: contrasting "resummed with 
Taylor [p]" approximants with some 5PN-accurate "resummed 
with Pade[p]" approximants. See text for definitions and ex- 
planations. 



wc discuss here only, besides the dominant even-parity 
quadrupole mode |/i22|, the first subdominant odd-parity 
mode |ft.2i|, as well as the dominant £ = A mode |/i44|. 

Fig. 6 focusses on |/i22|- The left panels, (a) and (c), 
display the standard Taylor-expanded |/i22| = 1 + cix + 
cx.^x^l"^ + cix^ -|- . . . . More precisely, panel (a) consid- 
ers the standard Taylor-expanded amplitudes up to 3PN 
accuracy included, while panel (c) displays the standard 
Taylor-expanded amplitudes from 3PN to 5.5PN accu- 
racy. By contrast, the right panels, (b) and (d), dis- 
play our new "resummed with Taylor[p]" approximants: 
panel (b) exhibits the IPN, 2PN and 3PN approximants, 
while panel (d) contrasts the 3PN, 4PN and 5PN approx- 
imants.^^ Consistently with previous studies [5, 6] (done 
at the level of the flux) there is evidently more scatter in 
the standard Taylor-expanded amplitudes than in the re- 
summed ones. In particular, note that the standard IPN- 
accurate Taylor approximant gives a grossly inaccurate 
representation of /i22 as soon as a; > 0.05 (building up 
to —40% at the LSO), while our new-resummed Ti[p22]- 
based waveform not only captures the qualitative behav- 
ior of the exact waveform, but also reproduces it quan- 
titatively within ^ 4% even at the LSO. On the other 
hand, for 3PN and higher accuracies the resummed wave- 



forms exhibit a very close agreement (within ±1 x 10 
with the exact one.^^ In previous work [1, 2], we had 
proposed to resum /122 by Padc (-P2) approximating 
/22 = (^22)^ instead of ^22- For completeness, we com- 
pare in Fig. 7 our previous best proposal to the cluster of 
our current best proposals (based on various Taylor and 
Pade approximants of Pti)- In first approximation, this 
figure shows a rather close agreement between all these 
approximants. In second approximation, one can note 
that some of our new approximants, namely P3 , P\ and 
Pj*, are closer to the exact numerical results. From the 
pragmatic point of view, our current best-bet approxi- 
mants are therefore our two new, pole-free, Pade approx- 
imants based on P3^{T5[p22]} and P\{T^\p22\\ ■ We have 
a slight preference for P^\T^\p2-2\\ which is the normal 
sub-diagonal Pade (admitting a simple continuous frac- 
tion representation) and which was close to the sub-LSO 
numerical results (see Fig 5). 

Fig. 8 exhibits the results for |/i2i|. We compare 
and contrast: (i) standard Taylor-expanded amplitudes 
(top panel), (ii) new-resummed amplitudes when factor- 
ing J (middle panel) and (iii) new-resummed amplitude 
when factoring TJoff (bottom panel). For brevity, the 
"standard Taylor-expanded" top panel exhibits only the 
integer-order PN-approximants. Note again, as in the 
case of \h2i\ discussed above, how the use of a stan- 
dard IPN-accurate Taylor-expanded waveform leads to 
a grossly inaccurate approximation to the exact result, 
building up to —22% at the LSO. By contrast, our new- 
resummed ri[/92i] "based approximant (middle panel) or, 
for that matter, the Ti [p|^]-based one (bottom panel), 
captures both qualitatively and quantitatively the cor- 
rect behavior of the exact waveform.^" Ultimately, for 
3PN and 4PN-accuracies, both resummed waveforms ex- 
hibit a very close agreement (within ~ 3 x 10"'^ for the J 
case) with the exact one. The standard Taylor-expanded 
ones are also close to the exact results, but visibly less 
close than our new approximants. 

Finally, Fig. 9 exhibits the results for the |/i44| wave- 
forms. Note that /144 and h^i are the last partial 
multipoles for which the analytical v = ^ result is 
known to 3PN accuracy. The comparisons between stan- 
dard Taylor-expanded and new resummed waveforms 
displayed in Fig. 9 leads to essentially the same con- 
clusions as above. In particular, the standard Taylor- 
expanded IPN accurate waveform is even more grossly 
inaccurate^^ than before, as the difference builds up to 



Fig. 1, between the exact total flux Fix) and the results ob- 
tained by using only "resummed with Taylor [p]" approximants, 
we will not discuss here the probable improvements that a further 
Padeing of the subdominant pirn's- might bring in. 
Note that because our tail factor (together with e"'*™ ) has conve- 
niently resummed all the half-integer powers of x, the left panels 
have to include half-integer PN-approximants, while the right 
panels have only integer-power approximants. 



As before, the fact that the resummed T3[p22] approximant is 
closer to the exact result than the T4[p22] and r5[p22] ones is 
probably coincidental. It is more important to note that panel 
(d) exhibits much less scatter than panel (c). 
The fact that the resummed Ti [p^il-based approximant is ex- 
tremely close to the exact result (see bottom panel) is proba- 
bly coincidental. We do not expect this coincidence to hold for 
higher-order partial waves. 

This is the consequence of the analytical fact noted above that 
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about -90% at the LSO (i.e. \hl^^{xLso)\ » 0.0793 
instead of \h^r''^\{xLso) = 0.8334)! Let us also em- 
phasize that, as we could have aheady pointed out for 
I/122I and for |/i2i|, the new resummed approximant based 
on the 2PN- accurate Taylor-expanded pem's is systemat- 
ically less good than the one based on the IPN-accurate 
Taylor-expanded ptm^s. This suggests that for wave- 
forms which are subdominant with respect to /i44 and 
/i42 (for which one does not know the 3PN expansion of 
the waveform) one will be better off, if one intends to use 
Taylor-expanded pgrnS, in employing only the IPN accu- 
rate pirnS. However, as we have shown in the £ = m = 2 
case (see Table II and Fig. 7), we expect that a suitable 
Padc resummation of the highest accuracy available re- 
sults will yield better agreement than simply using the 
Taylor IPN-accurate pem's. In this respect, let us recall 
that, as exhibited in Eqs. (45)- (46), the IPN corrections 
for all even and odd-parity multipoles are known. In the 
limit they are given by Eqs. (45) and (46); in 
the comparable-mass case the even-parity result is given 
by Eq. (41) while the odd-parity result is given in Ap- 
pendix A. 



IV. RESULTS FOR THE COMPARABLE MASS 
CASE (NOTABLY THE EQUAL MASS CASE, 

u = 0.25) 
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Let us continue to test our resummation procedure by 
considering the comparable mass case v =^ 0, and no- 
tably the equal-mass case, v = 1/4 = 0.25. In this case, 
we cannot rely on the knowledge of the "exact" multi- 
polar waveforms from comparable mass circular orbits. 
Indeed; though this problem can in principle be numer- 
ically investigated for binary black hole systems by con- 
sidering the helical Killing-vector approach (see Ref. [34] 
and references therein), there arc no presently available 
results where one goes beyond the conformally flat ap- 
proximation to Einstein equations. (But see Ref. [35] for 
the case of binary neutron star systems). For what con- 
cerns the available numerical results on coalescing black 
holes, previous work has shown that the deviations from 
the adiabatic-quasi-circular approximation were far from 
being negligible near the LSO, so that they cannot be 
directly compared to the circular waveforms investigated 
in this paper. We leave to future work a comparison 
between suitably non-quasi-circular corrected analytical 
waveforms and the results of numerical simulations of 
coalescing black holes. 

In the absence of exact waveforms to be compared to, 
we shall content ourselves here by investigating the inner 
consistency and robustness of our current best-bet ana- 



the IPN correction to the waveform is negative and grows linearly 
with £. We recall that this fact was one of our motivations for 
introducing the new quantities pim- 
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FIG. 8: Extreme-mass-ratio limit {y = Q): various represen- 
tations of the |/!.2i| waveform modulus. Top panel: standard 
PN-expansion. Middle panel: resummation factoring the an- 
gular momentum J . Bottom panel: resummation factoring 
the energy Heg. 



lytical approximants as suggested by the = results 
reported above. More precisely, we shall study the de- 
pendence of pim and the corresponding new resummed 
waveform |/i£„i| on the two EOB deformation parameters 
v and 05. 
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FIG. 9: Extreme- mass- ratio limit (y = 0): various representation of the |/i44| waveform modulus. Left panel: standard 
PN-expansion. Right panel: our new resummation. 



A. Mild dependence of on v For concrctcness, let us display here Pii{x; v) (given to- 

gether with our results in Appendix C) 

To motivate the study of the z^-dependence of pim let 
us start by having a close look at the general structure 
of Pirn, i.e. 

pim{x-v) = l+cf"(i/)a;+c^^'"(i^)x+cg^'"(log(x);i^)x3+. . . . 

(61) 



, , 2625i/2 - 58701/ -)- 1614 
=^ + 1320(3.-1) ^ 



1252563795.4 - 6733146000.3 - 313857376.^ + 2338945704. - 511573572 ^ 

317116800(1- 3.)2 ^ 
/ 16600939332793 12568 , , , A c, 

\ 1098809712000 ~ ^^"^^""^^S^^") j " ' ^^'^ 



We see on the example of P44 that the .-dependence of 
the coefficient c^'^™ (.) is not polynomial in ., but rather 
given by a rational fraction. The denominator of this ra- 
tional fraction in the case of P44 is proportional to some 
power of 1 — 3.. The denominator 1 — 3. decreases signif- 
icantly (from 1 to 1/4) as . increases from the extreme- 
mass-ratio case, . = 0, to the equal mass case, . = 1/4. 
From Eq. (7), for the general multipole this denom- 
inator would be proportional to a power of 

ai,n{i^) - ^2 + (-)'"Xi " X2 + (-)™Xi 

(63) 

This ratio is expressible as a polynomial in .. For in- 
stance, for P54, it would be c?54 = 1 — 5. -|- 5.^, which 



decreases from 1 down to 1/16 as . goes from to 1/4. 
More generally, dira{v) decreases, as . varies from to 
1/4, from 1 down to 1/2^+'"^ when 7r(m) = and to 
(£ + e — l)/2^+'^~^, when 7r(m) = 1. The presence of such 
"small denominators" raises the issue of a possible large 
increase of the coefficients cf^'^ (.) as . increases from to 
1/4. If that were true, this would undermine the applica- 
bility to the comparable mass case of the conclusions that 
we have drawn above from the . — > limit. Therefore, we 
have studied the .-dependence of the known coefficients 
c5^*'"(.) to check whether the presence of these "small 
denominators" might cause them to grow uncontrollably 
when . increases. In Table IV we list the fractional dif- 
ferences Ac^'!" (.) = c^'!" (i^)/<'" (0) - 1 at .= 1/4 (and 
at log(a;) = log(l/6) for the logarithms contained in the 
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TABLE IV: Analysis of the fractional variation 
Ac^*^" (ly) = c^'"- [u) /c?,''" (0) - 1 of the coefficients c^" (ly) 
in Eq. (61) for a selected sample of values of {£, m). 
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FIG. 10: Equal- mass case {u — 1/4): contrasting vari- 
ous methods for resumming the waveform modulus \h22\ for 
as = 0. Note the presence of a localized spurious pole in 
Pt{T5[p22]} atxKi 0.038. 



higher coefficients) for a sample of the pim's whose in- 
dependence is analytically known. The good news is that 
Table IV indicates that the fractional variation of the co- 
efficients c^^^ (v) when going from the extreme-mass ratio 
case to the equal-mass ratio one is typically of the order 
of 20%. 

This mild dependence of the coefficients c^^'™ (i/) on 
is the basis of the proposal [2] of improving the accuracy 
of known j/-dcpendent ptm^s by adding the v ^ limit 
of higher order PN-corrections ("hybridization"). [For 
instance, in the case of p22ix',h'), where the i^-dependent 
terms are known up to 3PN, we have added the 4PN and 
5PN ly ^ corrections]. Indeed, this procedure consists 
in using, for some higher corrections, the approximation^^ 
C'"(^)«<-(0). 



Note that our results on the mild ^-dependence of cj^*"* (u) show 
that, a contrario, a naive hybridization of the waveform of the 
type h"^^'"^ = h|°°""(in) + hY^''"(i' = 0) would probably be 
rather unreliable, especially for £ > 3, because it would not in- 
corporate the overall strong decrease approximately proportional 
to the "small denominator" d£„(j/), Eq. (63). 
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FIG. 11: Equal-mass case {y — 1/4): Same as Fig. 10, but 
focussing on only the "best" (3^'^PN-accurate) Fade approx- 
imants to the waveform. 



We have validated this approximate completion of 
known z^-dependent terms in the following way. In view 
of the results of Table IV, we have tested our procedure 
by modifying the 4PN coefficient for i/ = 0, C4^^(0), by 
multiplying it by the factor (1-1-0.8;^), in order to mimic a 
possible 20% increase of this coefficient when ly increases 
up to 1/4. We then found that such a modification of the 
4PN coefficient yielded a corresponding modification of 
T5[p22{x; ly)] equal to n[p22{x; /T^[p22{x; y)] - 

1.00038 when evaluated for v = 0.25 at a; = 1/6. Even 
at a; = 1/3 we find that such a modification yields 
Tz[p22{x-iy)]^°'''^°''m[p22{x;v)] = 1.013. In the p44 
case, where the i^-dependent corrections are known only 
up to 2PN accuracy, a similar modification of the 3PN 
term for = by a factor (1 -I- yields a corre- 

sponding fractional change of P44 between u — Q and 
V = 0.25 equal to 1.0099 at x = 1/6 and 1.079 at x = 1/3. 
These results confirm the reliability of the hybridization 
procedure adopted here, and give us an idea of the re- 
lated small uncertainty. For instance, for the dominant 
quadrupolar wave, we can anticipate that our hybridiza- 
tion procedure introduces an uncertainty in the waveform 
h22{x) cx {p22{x)f of order 8 x lO""^ at the LSO. This 
level of uncertainty is comparable to the fractional differ- 
ence between our best-bet quadrupolar amplitude based 
on P^[p22{x, V = 0)] and the exact result (see Table II). 

B. Mild sensitivity of \him\ to v 

In Fig. 5 we had put together, in the extreme-mass- 
ratio case, the predictions for |ft,22| made by all the 
higher-order approximants within our new resummation 
method. Let us now "deform" the results of Fig. 5 by 
turning on v and increasing it up to = 1/4. Fig. 10 
is the "i/ = 1/4-deformcd" version of Fig. 5. In con- 
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structing this figure we have used the value 05 = 
for the 4PN EOB parameter entering Eq. (11), and 
we have defined the EOB radial potential A{u) as be- 
ing P^[A^^y^°'-' {u; a^)]. The horizontal axis has been ex- 
tended up to the location of xi,so{o-5,'^) as predicted 
by the corresponding adiabatic EOB dynamics, namely 
a^LSo(0,l/4) = 0.2112. Some of the lessons we might 
draw from comparing the v = 1/4-deformed Fig 10 to 
its I' = counterpart. Fig. 5, arc the following: (i) 
apart from Pi{T^[p22]} (which still has a spurious pole) 
and our old -P2 1/22}, the relative stacking order of all 
the other approximants is maintained in the deforma- 
tion between v = and v = 1/4; (ii) our old prescrip- 
tion [1, 2] based on P2{/22}, which in the = case 
was clustered together with the other approximants (as 
well as with the exact curve), seems now to have drifted 
apart from the cluster of the other ones; (iii) indeed, 
all the new approximants are rather well clustered to- 
gether, with a dispersion which reaches only about 2% 
at XLso(0, 1/4) = 0.2112. 

One of the results of the 1^ = study above, partic- 
ularly in the dominant £ ^ m = 2 case, was to select, 
among the array of new approximants, a small sample 
of "best approximants" . This sample was made of the 
approximants based on P|{T'5[y022]} and P4 {Tr,[p22]} ■ In 
Fig. 11 we extracted from the previous figure the v = 1/4- 
deformed version of only these two "best approximants" . 
Remembering that in the 1/ ~ case these two curves 
were both extremely close (within 6 x 10~^ at x = 1/6!) 
to each other, as well as being very close to the cor- 
rect answer, we note that their v = 1/4-deformed ver- 
sions are still very close to each other (within 3 x 10"'^ 
at xlso(0, 1/4) = 0.2112 ). We therefore expect that 
this doublet of curves is a good indication of where the 
currently unknown (circular, adiabatic) correct v = 1/4 
curve might lie. 

C. Weak dependence of \him{x)\ on as. 

Finally, we study in Fig. 12 the sensitivity of our new 
resummed circular waveform to the 4PN EOB parameter 
05. This sensitivity comes from several sources. Both the 
source term sj^^ in Eq. (1) and the tail term Tgm depend 
on the EOB dynamical quantities H and J7. Therefore, 
when expressing the waveform as a function of the fre- 
quency parameter x, obtained by solving Eq. (13) above, 
the a5-dependent radial potential A{u) comes in at sev- 
eral different places. 

For concreteness, we shall study the "deformation" of 
our two best approximants when 05 increases from to 25 
(such a range is motivated by recent work [2, 10, 11, 18] ). 
The a5-dcformcd version of Fig. 11 is plotted as Fig. 12. 
This figure compares two doublets of curves: our two 
best Fade approximants {P^ {T5[p22]} , Pi{T5[p22]}) for 
05 = versus the same Fade approximants when 05 =25. 

The main thing to note is that the 05 deformation is 
continuous and monotonic. The displacement of each 
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FIG. 12: Equal-mass case (y = 1/4): EflFect of varying 05 (be- 
tween and 25) on the "best" Pade approximants displayed 
in Fig. 11. 



curve is only of order 2 x 10"^ at a; = 0.2112 = 
a^LSo(0, 1/4) (the horizontal axis of the figure has been 
extended up to a;Lso(25, 1/4) = 0.2236). In addition, the 
separation of the a5-deformed doublet of Pade curves is 
about the same as it was before deformation. 



V. CONCLUSIONS 

In this paper we have explored the properties of a 
new resummation method of post-Newtonian multipolar 
waveforms from circular nonspinning compact binaries. 
The two characteristic features of this method are: (i) 
the multiplicative decomposition of the (complex) /i^^ 
waveform into the product of several factors correspond- 
ing to various physical effects, and (ii) the replacement 
of the last (real) factor, /fm, in this decomposition, by 
its £-th root pe,nix) = {fim{x)Y^'^ . 

To test this resummation method we have first con- 
sidered the extreme-mass-ratio limit [v 0), for which 
"exact" results for the waveform can be obtained by nu- 
merical analysis of black hole perturbation theory. We 
first noted (see Fig. 3) that the new quantity that we 
introduced, pem{x), has a remarkably simple quasi-linear 
behavior as a function of the orbital frequency parameter 
X = (GMJl/c'^)^/'^. [In the odd-parity case, this quasi- 
linear behavior is especially pronounced when factoring 
out the angular momentum J7 from the wave amplitude. 
This leads us to consider p'^^^ as our "best-bet" default 
choice.] We related the simple properties of the func- 
tion pirnix) (including those concerning its dependence 
on {i,m)), to analytical results on the IPN corrections 
to multipolc moments. In this regard, we explicitly com- 
puted new expressions for the IPN source current multi- 
poles for arbitrary £ and in consequence the coefficient of 
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the IPN correction in the odd-parity waveform (and p£m)- 
The quasi-hnear behavior of the functions ptm{x) also 
means that 2PN and higher-order corrections to them 
are smaller than analogous corrections in usual quanti- 
ties, like the waveform. 

We have shown that, even if one uses only (without 
any further rcsummation) the successive Taylor approx- 
imants to pirm this defines a sequence of new resummed 
approximants to the waveform which "converges" to- 
wards the exact waveform much less erratically than the 
standard PN approximants. Moreover, for all the wave- 
forms for which 3PN corrections are known (at least when 
v 0), our results show that the new resummed wave- 
form nearly coincide with the exact results starting with 
the 3PN approximation (see Figs. 6, 8, 9). We have 
also shown that we can further improve the quality of 
our new approximants by suitably Padc-resumming the 
function p(m{x) before using it to construct the wave- 
form h'f^{x) (X {pim{x)Y ■ In particular, two Padc ap- 
proximants to P22, namely P|{T'5[/922]} and {Tr-,[p22]} , 
stand out as defining the most accurate representation of 
the exact waveform (see Fig. 7). 

We have finally explored the robustness of our approx- 
imants when considering a finite mass ratio. We have 
checked that the i^-depcndencc of the coefficients entering 
the Taylor expansion of the function pirn{x\v) is rather 
mild in spite of the presence of z/-dependent denomina- 
tors that decrease very significantly as v increases from 
to 1/4. This justifies the proposal of completing the 
known i/-dependent pem^s by adding the — > limit of 
higher order PN corrections. We have also shown that 
the relative stacking order of all the best approximants 
is maintained in the 'V-deformation" between v = Q and 
v — 1/4. In addition, our new approximants are rather 
well clustered together, with a dispersion which reaches 
only about 2% at the Last Stable Orbit. 

Let us finally note that we have compared in the four 
panels of Fig. 1 four different resummation approaches 
to the total (Newton-normalized) GW energy flux F{x) 
(for ^ 0): (a) the standard post-Newtonian (Tay- 
lor) expansion, (b) the Pade resummation advocated long 
ago [6] , (c) the improved ?;poie-tuned Pade resummation 
advocated in [2], and (d) our present new resummation 
method (using only Taylor expanded ptm&). The Wpoio- 
fiexcd technique, panel (c), is clearly superior to the re- 



sults of the first two techniques, panels (a) and (b). It 
has however the disadvantage that it needs to rely on 
some external knowledge (such as the exact value of the 
fiux at the LSO) to determine the optimal value of Upoic- 
On the other hand, our new resummation procedure not 
only stands out, among all other proposals, as yielding 
the best agreement with the exact flux (when ;/ = 0), but 
it has also the further advantage of being parameter-free. 
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APPENDIX A: RESULTS FOR ODD-PARITY 
(CURRENT) IPN-ACCURATE MULTIPOLES 

The IPN-accurate results for the source current (i.e., 
odd-parity) moments (of any multipolar order t) were ob- 
tained long ago in Eqs. (5.18) and/or (5.21) of Ref. [31]. 
Alternatively, we can use as starting point for the explicit 
determination of the IPN-accurate source current multi- 
pole moment Eqs. (4.3) and (4.4) in [36]. Recalling the 
notation 7 = GMjRc? and (consistently with Eq. (7)) 
using the notation 

h,(v) ^ Xi + {-fxi (Al) 
c,{v) ^ + (-)^A(-i (A2) 

,n^(yi,y2)^/^/^y^(yi,y2) (A3) 
oyi dyl 



1 ^ 

p=0 



at IPN accuracy the "compact" terms for the source cur- 
rent multipolc moment can be explicitly evaluated, in 
the circular orbit case, for a general value of £ (as in the 
t -\- m even case) . They read 



Tcompact 



STFl vMeabiiXaVbi x 



^(^-^^(-) + ^-^-"^(-)+(V 2l£;2!l2£V3) '--^-(-) 



+ —XL-3Vii_^Vi^_^ 



(£-l)(£-2)(^ + 4) 
2(£-h2)(2£-h3) 



(A5) 



r 



In addition to the above "compact terms" (generated by sor r'^'^), there exist three "non-compact" contributions 
compact-support terms in the effective stress-energy ten- 



24 



TABLE V: List of the bi{i') and ce{i') functions that appear in the text for some values of £. In the following formulas, we have 
introduced the notation X12 = Xi — X2 ~ sign(mi — m2)vT^-^- 



1 -X12 

2 l~2u 1 

3 -Xi2(l-!^) -X12 

4 1 - 4i/ + 2!^^ 1-3!^ 

5 -Xi2{l - + u^) -Xi2(l-2i.) 

6 l-ew + 9u^ -2u^ l-5u + rw'^ 

7 -Xi2{l-5u + 6u^ -u^) -Xi2{l-4u + 3u^) 

8 l-8iy + 20u^ - 16u-' + 21/" l-7u + 14,u^ - lu'-^ 

9 -Xi2{l-7u + l^u'^ -lOu^ + u^) ~Xi2(l-6i/ + 10;/2-4//3) 

10 l-lQu + S5u'^ - 50u-^ + 25!^" - 2u^ 1-91^ + 27u'^ - iQu^ + 9u* 

11 -Xi2(l - 9i/ + 28!/^ - 35i/^ + 15!/'' - i/^) -Xi2(l - 8;/ + 21i/2 - 20;/^ + 5!/-*) 



that make the IPN current moments more involved than 
the corresponding IPN mass moments. These noncom- 
pact contributions can be expressed in terms of the 
objects introduced in [31], so as to obtain. 



Tnoncompact 



GM 



2£+ 1 



L-lc 



L-la 



2{e + 2){2i + 3) dt 



( Y^' 



Icb 



^ 2 
(A6) 



More explicit expressions for these non-compact con- 
tributions can be provided for a general value of I by 
straightforward but slightly long computations. For cir- 
cular orbits one can check that the last term does not 



contribute and the final result for the other two terms can 
be simply re-expressed in terms of the polynomials hi{v) 
and ct{v), Eqs. (A1)-(A2), as for the compact terms. The 
final result (for circular orbits) is given by 



jnoncomp t:^ i\ t 



L-1 



C<;+3(z/) -I- 3&f+i(t/) 

2i 
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(A7) 



In the test-mass limit {v — > 0) this expression reduces to 



Tnoncomp 



STFl 2iyMjeabieXaVbX 



L-l(-l)' 



(A8) 



Thus, in the circular orbit case, the IPN-accurate current 
multipole for a general value of £ finally reads: 



Jl =STFl vMeabieXaVb<yX^ ^ 

lf£+l (^-l)(^ + 4) 



(£ + 2) (2^ + 3) 



ct+i{v) + 7 



^ , ^ 2£ + 3^ , , ^ ^-f 1 . , 
Q+i(i/) -I ;^^6£+i(i/) -I- 2v——bi^i\v)- 



2£ 



2£ "^'^ ' ' i 
i£-l){£-2)i£ + A) 



2{£ + 2){2£ + 3) 



(A9) 



Adapting the reasoning line of Rcf. [12]. recalling the 
additional velocity dependence of the current moments 

that leads to {v/cY^^ and noting that 7 = .t to this ^(^i-^ ( / 

order of accuracy, one can finally show that, for circular ^imi-'^^ i') ~ 1 — x<. {£ + 1) (1 

------ -- ni), - 



3/ ^ 2£ 



orbits, the IPN-accurate odd-parity h^J& read 3^ -h 3 hi+i{v) £ + 1 h_i{v) 

21 Q+iH"^'' £ c,+iM ^^^"^^ 
1 £-Hq+3M m^(£ + 4) c,+z{v) \ ^ 
2 £ Q+i(:/) + 2(^ + 2)(2£-f 3)c,+i(z.)J ^ ^' 
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where we have not simplified on purpose in order to allow 
the reader to explicitly track the origin of each single con- 
tribution. In the extreme-mass-ratio limit. M = mi ^ 
fj. = m2 {i^ = m/M 0), one has q(0) = 6^(0) = (-1)^ 
and so this equation simply reduces to 



1 



m2(£ + 4) 
2{e + 2)(2e + S) 



(All) 



When computing the amplitude fijn{x;v) (where 



(1,./) 



j is factorizcd), an additional contribution of 



-(3/2 + iy/6)x (see Eq. (18)) comes in, so that the IPN- 
accurate fj/,y, {x; i')'s read 



//™(x;.)=i-J(^ + l)(i-0+^ + ^ + ^ 



2€ Ci+i{v) 



2v- 



2 I c,+i(i/) ^ 2(£ + 2)(2£ + 3)Q+i(;.) 



+ Oix^). 



(A12) 



r 



In the test-mass limit, this equation becomes 



//„(x;0) = l-x £+1-- + 



+0{x^). 



These results lead to the IPN-accurate p 
and (49), that we have used in this paper 



(A13) 
/„;s, Eqs. (46) 



For completeness, we conclude this Appendix by quot- 
ing the I'-dependent, IPN-accurate fim^s for £ + m even 
and a ready-reckoner of the bi{i/) and Ci{i') functions for 
^-values relevant for this work. Table V. From Eq. (C5) 
of Ref . [12]. the general expression of hfj^ at IPN reads 



2 ci{v) ceiiy) ci(v) 2(£+l)(2£ + 3) 



+ 0{x?). 



(A14) 



From this expression, the even-parity fem^s follow as 



fe^{x;,.)^l-xil{l-^) 



1 3 ce+2{iy) ^ bt{v) ^ ct+2{v) m^(^ + 9) 

2 2 ci{v) 



Ci{v) ci{v) 2(^ + l)(2^ + 3) 



+ 0(x2). 



(A15) 



r 



which reduces to Eq. (43) in the test-mass limit. 



APPENDIX B: EXPLICIT FORM OF THE /f„'S 
WITH HIGHER PN ACCURACY 



multipoles up to = 5 for even-parity modes {£+m even) 
and ^ = 4 for odd-parity modes {£ + m odd) . We consider 
separately the even-parity /fm's and the odd-parity f^^^s, 
and fP^s. 

1. Even-parity fim's 



In this Appendix we complete the information given 
in the text by explicitly listing the fi,nS that are known 
at an accuracy higher than IPN. This means considering 



The even-parity firnS (with v and v ^ contri- 
butions) are given by 
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^ , , ^ 1 - 67451/ - 4288 2 
/22 (a;; = 1 + — (551. - 86)x + — 



/ 1146351/3 227875i^2 4^ 34625i/ 856 , , , , 21428357\ , 
1 u eulerlognicc H x 

V 99792 33264 96 3696 105 62 727650 / 

/36808 , , , , 5391582359\ , /458816 , , , ^ 93684531406\ . 
eulcrlogo(x) a; + eulerlogofx) + Ofa;^), (Bl) 

V 2205 ' 198648450 J \ 19845 ^ 893918025 / ^ ^ v ; 

7\ /8871/2 3401J, 443\ /147471561 78 , , , A , 



h,[x- .) = 1 + ( 2. - - j X + - W - 440 j " + 1^80280^ " yenlcrlog3(x) J . 

39 eulerlog3(x) - ^|^) + O(x^), (B2) 

, , , { 2u li\ ( 247t/2 i273\ , 

^^^("'") = ^+(,-y-yj"+(,-^-T9^ + wj" 

/400427563 26 . n\ 3 /'ISg 12064573043\ 4 5 



V 75675600 ^x'^^^^^'-^^^-^ ) + ^^-l-log,(.) - j^^^^ ) + Oix^), (B3) 



r 



J 



, , , 26251.2 - 58701/ + 1614 23740185i/3 - 106831480i/2 + 50799672i/ - 4536144 2 

/44 UC; J/ = 1 H ^ ^ a; H ^ x 

^ ' ' 330(3i/-l) 1801800(3z/- 1) 

1132251120 ^ , , , 5992751383\ 3 , 4, 

culcrlog.(x) ] x-^ + 0(x ), (B4) 

156080925 ^ 156080925 ^ ^ ^' V ; 

285t/2-3530t/+ 1146 2707215:/^ + 28154560;/^ - 26861688;/ + 5538096 2 
/42(a;; z/) - 1 + 330(3^, _ 1) ^ 1801800(3i/ - 1) 

1132251120 , , , , 5180369659\ 3 4, 

eulerlogofx) x"^ + O(x^), (B5) 

312161850 312161850 J ^ ^ ' 



78(2t/ — 1) 6552 

(B6) 

^ , , ^ 176^^- 850^ + 375 69359 2 3^ 
■'^'^^ ' ' 78(2i/-l) 10920 ^ 

(B7) 

fzi{x- v) = 1+ 8^2 — 626;/ + 319 ^ ^ ^8859 ^2 ^ 0(2;3) L^t us focus now on the odd-parity case and list the 

' 78(2;/ — 1) 4680 ' f'^^ in which the Newton-normahzed angular momentum 

(B8) 

j has been factorized as an effective source. We have 
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, , f2iv 59\ /85i^2 269i/ 5\ , /88404893 214 , , , ^ 

fni (x: I') = 1 + \x + \ a; + eulcrloffn (a;) x'^ 

■'^^^ ' ' V 42 28/ V 252 126 9/ ^11642400 105 7 

/6313 , , , , 33998136553\ , 

+ - 4237833600 ) ^ + )' ^^^^ 

J 320;^2_ ;lii5j^_j^32s 39544^/^ _ 2537681^2 + 117215:/ - 20496 2 

/32(.t; - 1 + 90(3^ ^ + 11880(3:.- 1) 

/ 110842222 104 , , , A o 4, , , 

+ (^72972^ - ^-1-I°g2(-) j -3 + O(x^), (BIO) 

J, , (l60:/2 _ 547^ ^ 222) 225543 , , , , 

//3(^; ^) = 1 + ^ 44(2.-1) + 1004^^^ + 

//,(x; .) ^ 1 + 1385. + 602) 760181 

■^^^^ ' 132(2.- 1) 120120 ^ > ^ ' 



3. Odd-parity /i^'s have 

We finally list the odd-parity f^^^ in which the effective 
energy Hcff has been factorized as an effective source. We 



^H, ^ /'5. 3 \ /79.2 485. 97 \ , /70479293 214 , , , x\ s 

107 , , , , 9301790917\ , 

euler\og,{x) + ,,,,,,,,,, a;^ + 0{x'), (B13) 



490 t,iw ■ i4i2611200y 

^rr, , 365.2 - 590. + 148 16023.^ - 93976.2 + 612. + 6192 

fi;(a;: .) = IH ^ ^ x H ^ x'' 

■'^^^ ' ' 90(3.- 1) 3960(3.- 1) 

/96051082 104 , ^ , ^^ 3 4. .-o. .^ 

+ ll72972^ - ) + ^(-')' (B14) 

524.2 - 1135. + 402 i667 
/43(^,^)-l+ 132(2.-1) ^ 3640-'^+^^''^' 
, 332.2 - 879. + 338 145021 „ 

^) - 1 + 132(2.-1) - + 120120^ + '^(^^^ 



APPENDIX C: COMPLETE EXPRESSIONS OF 

THE pcm'S FOR 2 < ^ < 8 

We finally list the "hybridized" expressions of all the 
even- and odd-parity pirnS obtained from the corre- 
sponding /frn's [with the proviso explained above that 
the ^-th power of the "hybridized" pim presented here 



would generate some specific .-dependent higher-order 
coefficients c{5" (.) which differ from the c^*" (. = 0) 
listed in, e.g. the equations of Appendix B. In the odd- 
parity case, we only list the j/-normalized p/^'s obtained 
from the /^'^'s. For completeness and future reference we 
present the pim^s explicitly up to £ = 8 included. 



'bbv 43 Y / 195831/2 330251/ 20555 \ 
^"ST ~ 42 j ^ ^ V 42336 21168 ~ 10584 J 



/ 10620745^/3 6292061t/2 4l7r2t/ 48993925i/ 428 , , , , 1556919113\ , 

V 39118464 3259872 192 9779616 105 ' 122245200 / 

/9202 , , , , 387216563023\ , /439877 , , , ^ 16094530514677\ . 

culcrloKn ( X' ) \ X + \ eulcrlogo a;) x + Oil 

V2205 ' 160190110080/ V 55566 ' 533967033600/ ^ 



1+1^ 59\ ^ ^ / 617i/2 10993;/ 47009 \ ^ 



84 56 / ^ V 4704 14112 56448 / 



„2 



/ 7613184941 107 , , . ^^ ^ /6313 , , , , 1168617463883\ 4 

I I \ ^ ^ I — I'-k^'-'ll ^nl ,-1-1- \ /~i i-i- I rvt \ l/-y^^l/ l(/-^^'-'^ 



, culcrlogi (x) \ X + \ culcrloff, (a;) x + Oix ), 

V2607897600 105 ' J V5880 ' 911303737344 / ^ ^' 

, f2v 7\ /149^2 i^Qi^ 6719\ , /3203101567 26 , , . ^^ 

1 + ] x + [ a; + eulerlogofx x"* 

V3 6/ V 330 990 3960/ V 227026800 7 7 

13 , , , , 57566572157\ 4 5, 
-eulerlog3(x)- 8^62153600 j ^ + ^^^ ^^ 

320;/2 - 1115^ + 328 3085640^'' - 20338960^/^ - 47256051/^ + 8050045;/ - 1444528 2 
^ ^ 270(3;/- 1) ^ ^ 1603800(1 - 3;/)2 ^ 

/5849948554 104 , , , A 3 4^ 
(^4035532^ -^^^^^^-l^S^^^O^ 

2;/ 13\ / 829;/2 1685;/ 101 \ , / 11706720301 26 , , , ,\ 
Y-m)^+[-l782-l782 +n28)''+ [ 6129723600 ~ ^^"^^^-l^Si^^) j ^ 
/169 , , , , 2606097992581\ 4 5. 

2625;/2 - 5870;/+ 1614 
^ 1320(3;/- 1) 
1252563795;/^ - 6733146000;/3 - 313857376;/2 + 2338945704;/ - 511573572 

317116800(1- 3;/)2 ' 

/ 16600939332793 12568 , , , C\ 

eulcrloff4(a;) x"^ + O(x^), 

V 1098809712000 3465 V ' v 

, 160;/2 - 547;/ + 222 6894273 



.2 



a;- +0(a; 



176(2;/ - 1) 7047040 



285;/2 - 3530;/+ 1146 
^ 1320(3;/- 1) ^ 
-379526805;/-* - 3047981160;/^ + 1204388696;/^ + 295834536;/- 114859044 



317116800(1- 3;/)2 

/848238724511 3142 , ^ / n\ 3 , .0/ 4n 
eulerlog2(a;) x"* + 0{x^), 



x^ + O(x^) 



V 219761942400 3465 
, 288;/2 _ 1385^ _^ 602 7775491 



528(2;/ -1) 21141120 
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, , 512(^2 _ 2298z/ + 487 3353747 , ,r. . 

^-(^^ ^) = ^ + 390(2.-1) - - 2129400^ + ^^^"^ 

, , 33320i^3 _ ;l27610!^2 _^ 960191^ - 17448 2^ 
^) = ^ + 13650(5.2- 5. + 1) + )' 

, , , 176.2 _ 850. + 375 410833 2 3^ 

^) = ^ + 390(2.- 1) - - 709800^ + ^^l^) 

, 21980.3 - I04930i^2 _,_ g4679j^ _ 15828 2^ ^ 

^) - ^ + 13650(5.2- 5. + 1) ^ + )' (^1^) 

P5i(x; .) ^ 1 + ^-^-626. + 319 _ ^ 

' 390(2.- 1) 304200 ^ ^ ^ ^ 

, , , 273.3- 861.2 + 602.- 106 2^ 
P6a(x; -) = 1 + 84 (5.2- 5. + 1) ^ + ^^^^^ 

J/ ^ 1 220.3- 910.2 + 838.- 185 ^^,3, 
= 1 + 144 (3.2- 4. + 1) ^ + ^(-^ 

, , 133.3- 581.2 +462.- 86 3^ 
Pe4(x; .) = 1 + 84 (5.2- 5. + 1) ^ + )' 

J/ ^ 1 156.3 - 750.2 + 742.- 169 ,^,2, ,^,on 
Pa^3(-; = 1 + 144 (3.2- 4. + 1) ^ + '^(^ ) 

^ ^ ^ 49.3 - 413.2 + 378.- 74 2n .r^mN 

Pa2(x; .) = 1 + 84 (5,2 _ 5, + 1) ^ + (C19) 

7. N , 124.3 - 670.2 + 694.- 161 2n 
Pa^.(x; .) = 1 + 144 (3,. _ 4,^1) + (C20) 

, , ^ 1380.3 - 4963.2 + 4246.- 906 2^ .^^.s 

Prr(x; -) = 1 + 714 (3.2- 4. + 1) ^ + )' ^^^1) 

, , , , 6104.4 - 29351.3 + 37828.2 - 16185. + 2144 

^) - ^ + 1666 (7.3 - 14.2 + 7.- 1) ^ + )' (^''^ 

, , , 804.3 - 3523.2 + 3382.- 762 2, 
P75(x; .) = 1 + 714 (3,. _ 4,^1) + (C23) 

T, , , 41076.^^- 217959.3 + 298872.2 - 131805. + 17756 2n /^^.n 

^) - 1 + 14994(7.3 - 14.2 + 7.- 1) ^ + ^^^^^ 
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420i/3 - 2563J/2 + 28061/ - 666 , ^ 



P73(x; .) = 1 + (3,. ^4,^1) + 0[.% (C25) 

, , 32760(^4 - 190239!^3 + 273924(^2 _ ;L23489i/ + 16832 

^) = ^ + 14994(7.3 - 14.^ + 7.-1) ^ + ^^^^^ 

, , , 228.3- 2083.2 + 2518.-618 2^ 

Pn(x; .) = 1 + (3,. _ 4,^1) + (C27) 

, , 12243.4 - 53445.3 + 64659.2 - 26778. + 3482 /^^on 

^) - ^ + 2736(7.3 - 14.2 + 7.-1) ^ + '^^^ ^^^^^ 



, , 38920.-4 - 207550.3 + 309498.2 - 154099. + 23478 2^ 

Pc7(x: .) = IH ^ — = ^ ^ X + 0(2:2), /^gg) 

'^^^^ ' ' 18240(4.3 - 10.2 + 6.- 1) V V ; 



, , , 2653.4 - 13055.3 + 17269.2 - 7498. + 1002 2x /^o^x 

P86(a;;.) = l + , — , , ^ -r x + C'(x2), (C30) 



P8^5(x;.) = l + 57T7F7T^ ,n 2 , x + 0(x2), (C31) 





912(7.3-14.2 + 7.-1) 


X 


6056.4 


- 34598.3 + 54642.2 - 28055. 


+ 4350 




3648 (4.3 - 10.2 + 6.- 1) 




4899.4 


- 28965.3 + 42627.2 - 19434. 


+ 2666 



^-(^^ ^) = ^ + 2736 (7.3 - 14.2 + 7.- 1) ' + )' ^^^^^ 

J, , , 24520.4 - 149950.3 + 249018.2 - 131059. + 20598 
^-(-' = 1 + 18240 (4.3 - 10.2 + 6.- 1) ^ + 



, , , 3063.4 - 22845.3 + 37119.2 - 17598. + 2462 

Ps2ix;i^) = l + 3 7T x + Ox2, C34 

2736 (7.-^ — 14.^ + 7. — 1) 



T , , 21640.4 - 138430.3 + 236922.2 - 126451. + 20022 

^) - ^ + 18240(4.3 - 10.2 + 6.- 1) ^ + 
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